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Abstract 

Membranes in M-theory are expected to interact via splitting and joining processes. 
We study these effects in the pp-wave matrix model, in which they are associated with 
transitions between states in sectors built on vacua with different numbers of membranes. 
Transition amplitudes between such states receive contributions from BPS instanton con¬ 
figurations interpolating between the different vacua. Various properties of the moduli 
space of BPS instantons are known, but there are very few known examples of explicit 
solutions. We present a new approach to the construction of instanton solutions interpo¬ 
lating between states containing arbitrary numbers of membranes, based on a continuum 
approximation valid for matrices of large size. The proposed scheme uses functions on 
a two-dimensional space to approximate matrices and it relies on the same ideas behind 
the matrix regularisation of membrane degrees of freedom in M-theory. We show that 
the BPS instanton equations have a continuum counterpart which can be mapped to 
the three-dimensional Laplace equation through a sequence of changes of variables. A 
description of configurations corresponding to membrane splitting/joining processes can 
be given in terms of solutions to the Laplace equation in a three-dimensional analog of 
a Riemann surface, consisting of multiple copies of connected via a generalisation of 
branch cuts. We discuss various general features of our proposal and we also present 
explicit analytic solutions. 
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1 Introduction 


The best candidate for a microscopic description of M-theory is the matrix model originally 
proposed in as a regularisation of the supermembrane theory and subsequently conjectured 
in to describe the full dynamics of the theory in light-front quantisation when the size of 
the matrices is sent to infinity. The matrix model provides a formulation of M-theory which 
is capable of describing multi-membrane configurations, arising as block-diagonal matrices. 
Membranes in M-theory should interact via splitting and joining processes and therefore the 
matrix model should capture such effects. However, no concrete proposal for the description of 
these splitting/joining interactions in the matrix model has been formulated. More generally 
there is no quantitative prescription for the study of this type of membrane interactions, 
which would make it possible to evaluate the associated transition amplitudes. This makes 
it difficult to test any results for these amplitudes that may be obtained from the matrix 
model. This is of course a general difficulty with all predictions of the matrix model, which 
so far have been mostly tested at the level of the low energy supergravity approximation. 

The AdS/CFT correspondence - and specihcally the duality proposed in , relating M- 
theory in AdS 4 x to an AA = 6 Chern-Simons theory - provides in principle new means 

of testing the matrix model predictions by comparing them to a dual CFT. In the context of 
the duality of it was recently shown in that the matrix model description of M-theory in 
AdS4x57Zfc can be quantitatively compared to the dual gauge theory without relying on the 
supergravity approximation or compactification to type IIA string theory in ten dimensions. 
The crucial observation of was that, focussing on large angular momentum states in 
M-theory and the dual CFT sector involving monopole operators, natural approximation 
schemes arise on the two sides of the duality, so that a systematic, quantitative comparison 
is possible. 

On the gravity side, M-theory states with large angular momentum, J, along a great 
circle in 5^ can be studied using the pp-wave approximation. The associated matrix model 
uses J X J matrices and its action was constructed in [^. We will use basic properties of 
the model which were further studied in [^. Multi-membrane states in the pp-wave matrix 
model consist of concentric membranes and their fluctuations. More precisely the vacua of the 
theory consist of spherical membranes which extend in AdS 4 directions and are point-like in 
5^. They are classified by a set of integers, Jj, i = 1,..., n, corresponding to a partition of the 
total angular momentum among n membranes. The fluctuations of the spherical membranes 
described by the pp-wave matrix model are associated, in the dual gauge theory, with certain 
monopole operators. The latter are characterised by their integer GNO charges, which are in 
one-to-one correspondence with the angular momenta, J^, of the membranes. The sector of 
monopole operators with large GNO charges can be reliably studied using a weakly-coupled 
effective low-energy approximation. 

The AdS 4 /GFT 3 duality in this M-theoretic regime relates correlation functions of 
monopole operators in the ABJM theory to processes on the gravity side in which the dual 
states interact in the bulk and propagate to the boundary. The simplest such process involves 
a single membrane splitting into two - with the associated three states propagating to the 
boundary - and it corresponds to a three-point correlation function of monopole operators. 
Similar processes have been studied in the case of the AdS 5 /GFT 4 correspondence, in which 
three-point functions of gauge-invariant operators are related to cubic interaction vertices in 
the bulk associated with the splitting of strings in the pp-wave approximation. A method 
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to compare the cubic vertices in string theory to the corresponding conformal three-point 
functions in the CFT was discussed in , following a spacetime interpretation first proposed 
in Studies using other approaches can be found in j9 -14 . Similar aspects have been 


extensively studied beyond the pp-wave approximation taking advantage of the integrability 
appearing in the analysis of the system. For recent interesting developments and up-to-date 


references, see 15,16 


The three-point correlation functions of monopole operators, in particular those with large 
charge J, which are relevant in the present case, are well defined in the ABJM theory and one 
expects that they should be computable within the CFT framework. Although in this paper 
we focus on the gravity side of the correspondence, it is important that in principle the results 
of the matrix model calculation can be verified by independent means. Successfully comparing 
the physical transition amplitudes associated with a membrane splitting or joining process 
to the gauge theory results for three-point functions would provide such an independent 
test. The agreement would not only represent a highly non-trivial test of the AdS 4 /CFT 3 
correspondence, but it would also provide strong evidence indicating that the matrix model 
captures the aspects of the dynamics of membranes related to splitting/joining interactions. 

In the context of the pp-wave matrix model, processes of joining or splitting of membranes 
correspond to transitions between states built on vacua with different numbers of membranes. 
A semi-classical description of these transitions can be given in terms of tunnelling amplitudes 
associated with instanton configurations which interpolate between states containing different 
numbers of membranes. In this paper we initiate a study of these interpolating configurations. 
The relevant instantons are BPS solutions to the Euclidean equations of motion of the pp- 
wave matrix model 17 , with boundary conditions specifying the single or multiple membrane 
states between which they interpolate. The equations obeyed by these BPS instantons, with 
the same boundary conditions, were studied in a different context in [^, where they were 
found to arise as the equations describing domain walls interpolating between distinct isolated 
vacua of the so-called N = 1* SYM theory. A number of properties of these instanton 
equations and the associated moduli spaces were studied in 18 . An important ingredient 


in the analysis presented in [18| is the observation that the relevant instanton equations can 
be mapped to the Nahm equations which arise in the construction of BPS monopoles, 
although with different boundary conditions. 

The BPS instantons of the pp-wave matrix model are saddle points representing the 
dominant contribution to the tunnelling processes. Physical transition amplitudes associated 
with the splitting/joining processes can be computed using a semi-classical approximation 
around the instanton solutions. This involves evaluating the determinants arising from the 
non-zero mode fluctuations as well as the integration over the collective coordinates associated 
with zero modes. Although there may be alternative approaches to the computation of 
the splitting/joining transition amplitudes, in order to carry out the standard semi-classical 
calculations it is necessary to obtain the explicit form of the classical instanton solutions. 
However, although various general features of the instanton equations - including conditions 
for the existence of solutions and properties of their moduli spaces - were studied in [18| , 
little is known about explicit solutions. In particular, no solution is known for the most 
elementary process we are interested in, i.e. the splitting of a single membrane into two. 

In this paper we show that an efficient approach to the construction of solutions to the 
instanton equations as well as to the study of their moduli spaces can be developed using 
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a continuum approximation. Such an approximation scheme is valid in the case of large 
matrices and, therefore, it is applicable in the context of the M-theoretic regime of the 
AdS 4 /CFT 3 duality considered in [^. In this approximation the original instanton equations 
can be mapped to a continuum version of the Nahm equations, which in turn can be shown 
to be equivalent to the three-dimensional Laplace equation, using results in 20-22 . In this 


formulation equipotential surfaces for the solution to the Laplace equation directly represent 
the profile of the membranes at different times. Using this approach we present explicit 
analytic solutions describing the splitting of membranes. 

A remarkable feature of our construction is that, in order to obtain a solution representing 
the splitting of one membrane into two, we are led to consider the Laplace equation not in the 
ordinary three-dimensional Euclidean space, but rather in a three-dimensional generalisation 
of a two-sheeted Riemann surface, which we refer to as a Riemann space following [23| . 
The emergence of Riemann spaces, which can be naturally motivated in the context of the 
construction of our solutions, is quite intriguing. The possibility that in the description 
of membrane interactions Riemann spaces may play a central role, similar to that played 
by Riemann surfaces in string perturbation theory, is very interesting and deserves further 
investigation. 

This paper is organised as follows. In section we describe the instanton equations 
and we recall some general properties of their solutions and of the associated moduli space 
which were obtained in |18| . In section we discuss the continuum approximation scheme 
used in our analysis and the mapping of the instanton equations to the three-dimensional 
Laplace equation. In section we present explicit solutions to the Laplace equation which 
are relevant for the pp-wave matrix model and we describe the corresponding solutions to the 
BPS instanton equation within our continuum approximation. The most interesting solution 
- which utilises a two-sheeted Riemann space to describe the splitting of a single membrane 
into two - is considered in section 14.21 We conclude with a discussion of our results in 
section Technical details are discussed in various appendices. 


2 BPS instantons 

In our study of configurations describing membrane splitting/joining processes we focus on the 
pp-wave matrix model . As discussed in this model provides a good approximation to the 
dynamics of M-theory in AdS 4 x in a sector containing membranes with large angular 

momentum, J. M-theory in this background is dual to the M = 6 Chern-Simons theory 
with gauge group U(A^)x \J{N) and level k constructed in [^. The region of applicability 
of the pp-wave approximation is defined by the conditions 1 <C J <C Moreover for 

J S> the model is weakly coupled. 

In this section we review the BPS instanton equations of the pp-wave matrix model and 
recall some of their properties. These equations were first studied in a different context in [18| , 
where they were shown to describe supersymmetric domain walls interpolating between vacua 
of the M = 1* SYM theory. Their relevance in the pp-wave matrix model was pointed out 
in [^, where they were first shown to describe instanton configurations interpolating between 
different vacua. 
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2.1 Instanton equations 

The Euclidean action describing the M-theory side of the AdS 4 /CFT 3 duality in the large J 
sector takes the form [4||^ 


Se = / dt tr< — ( - 

k 
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DX^ 

~15i 

2 ki 
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+ {2ETf— {{i[X^,X^]f + 2{i\X^,Y%‘^) 
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where T = l/[(27r)^/p] is the membrane tension, R is the radius of (and twice the radius 
of AdS 4 ) and 7 " (a = 1, 2,..., 9) are SO(9) gamma matrices, with 7 ^^^ = 7 ^ 7 ^ 7 ^. Following 
the notation in j^, we have denoted by T* (i = 1, 2, 3) and X"^ [m = 4, 5,..., 9) matrices 
originating from the membrane coordinates in AdS 4 and respectively. T is a 16 component 
matrix valued spinor. The y*’s, X^'s and the components of T are K x K matrices, with 


K = J/k. The covariant derivative in (2.1) is dehned by 


DX dA 


( 2 . 2 ) 


where Aq is the gauge potential associated with the invariance of the model under time 
dependent unitary transformations. In the following we will choose the gauge Aq = 0, which 
is compatible with the boundary conditions relevant for the tunnelling processes we are 
interested in. 


The matrix model dehned by the action (2.1) is a special case of a family of pp-wave 


matrix models, which are usually parameterised by a mass scale jjL. In the following we use 


the model (2.1), which is the one relevant for the large J sector of the AdS/CFT duality 


of and corresponds to a specihc choice of fi. To simplify notation and equations, we also 
focus on the k = 1 case. The generalisation is straightforward and in appendix [F| we present 
the main formulae for the case of a general pp-wave background. Using the equations in the 
appendix it is easy to recover the results for arbitrary k. 


The potential term in the action (2.1) is written as a sum of squares and thus it is 


manifestly non-negative. This makes it easy to identify all the vacua of the model, i.e. the 


zero-energy minima of (2.1). They correspond to the configurations 


{2ttT)R^ 


L\ with [L\ L^] = 


Ao™ = 0 , 


4'n = 0 


(2.3) 

(2.4) 


and therefore they are classihed by (generally reducible) J dimensional representations of 
SU(2). The solution in which the y*’s are proportional to the J dimensional irreducible 
representation corresponds to a single spherical membrane. It extends in the AdS 4 directions, 
with radius 

(2.5) 


r = 


{2'kT)R^ {2t:T)R? 
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and carries momentum J/R along a great circle of 5^. Generic vacua are given by block 
diagonal F* matrices corresponding to reducible SU(2) representations. The blocks have size 
Ji, i = 1,2,... ,n, with ^2- Ji = J and they represent n concentric membranes of radii 


{2-kT)R^ ’ 


( 2 . 6 ) 


with angular momenta Ji along the same great circle in 5^. 

We will denote the irreducible J dimensional SU(2) representation by J and the re¬ 
ducible representation which is the direct sum of irreducible representations of dimension 
Ji, J 2 ,. ■ ■, Jn by ^©^©•••© Jn - 

We are interested in the tunnelling processes corresponding to classical solutions inter¬ 
polating between a vacuum associated with the SU(2) representation in the inhnite 

past {t = — 00 ) and another vacuum associated with at t = + 00 . These processes are 

governed by the path integral with boundary conditions 


W(-oo) 

r(+oo) 


(27rr)i?2^(-°°) 

2 • 

{2t^T)R^ ^^(+ 00 )^ 


(2.7) 


In (2.7) we have chosen the generators to correspond to the standard embedding of 

SU(2) consisting of block diagonal J x J matrices. U denotes an arbitrary unitary matrix, 
which we need to include to take into account the gauge choice, ^0 = 0. 

In order to obtain the equations obeyed by the classical interpolating configurations we 
rewrite the bosonic part of the Euclidean action (2.1) in the form of a sum of squares plus a 
boundary term, 


Si; =^/dnr 


dX™\" 1 (27rr)27?2 


dt 


+ + 




+ (^^yW* + z(27rr)|eb^y*[y^ y'^]^ 

( 2 . 8 ) 


+ fermions 


From (2.8) it follows that, on configurations interpolating between vacua at t = ± 00 , the 
bosonic part of the Euclidean action obeys the bound 

1 


Se>^ tr( ^y*y* + i^^^eb''y*[y^y'=' 
R^ 6 


-1-00 


(2.9) 


Using the explicit form of the vacuum configurations at t = ±00 the bound (2.9) becomes 

, 2 

( 2 . 10 ) 
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T—iy[y] 


+00 


1 




^3i?2 (^(27rT)i?2 j ^’^(^(+oo)^Uoo) ^(-oo)^(-oo)^ 

where following we dehned the functional 

iy[y] = tr + i(27rr)^eb^y*[y^y^]^ . (2.11) 
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The BPS (anti-)instantons are the configurations m\ which saturate the bound ( |2.10[ ) 
and therefore satisfy the equations 


Ayi 2 7? 

— ± -y* ±^(27rT)-e*^■"[y^y^] = 0, 


( 2 . 12 ) 


with y™ = 0 and T = 0. In this paper, we refer to solutions of (2.12) in which one chooses 


the upper or lower signs as instantons or anti-instantons, respectively. We note that the 


Gauss law constraints are satisfied as a consequence of (2.12) 

dy* 


y* 


dt 


= T^( 27 r^):^e*^■^[y^ [y^y'^]] = o. 


(2.13) 


Since the (bosonic part of) the Euclidean action (2.1) is non-negative, so should be the ex¬ 


pressions on the right hand side of (2.9) and (2.10), when the bounds are saturated. Therefore 


iy[y] should decrease from t = —oo to t = -|-oo in the case of instantons (upper signs in (2.9) 


(2.12)) and increase in the case anti-instantons. This can be seen explicitly noticing that the 


instanton equations imply 




dW[Y] 


9y* 


(2.14) 


where the signs are correlated with those in (2.12). 


In the following we focus on instanton configurations and all the equations that we present 


correspond to the choice of upper signs in (2.12). Anti-instanton solutions can be simply 


obtained from the corresponding instanton by changing t to —t. 

A detailed discussion of the general conditions for the existence of solutions to the instan¬ 
ton equations (2.12) can be found in [l^. We will revisit these issues in section 4.3.2 Here 


we only recall a necessary condition. For a BPS solution connecting two vacua to exist, the 
vacuum with the larger value of iy[y] should correspond to a representation which does not 
contain more irreducible blocks than the one with the smaller value of iy[y]. This allows us 
to identify instanton conhgurations as corresponding to membrane splitting processes, while 
anti-instantons describe joining processes. 


Evaluating iy[y] on a vacuum conhguration (2.3) containing m membranes with angular 
momenta Jj, z = 1,..., m, one gets 






= . (2.15) 


Therefore on a generic instanton interpolating between the representations • • • © Jm 

and J 2 ® ■ ■ ■ ® J'n (with = J, n > m) the Euclidean action is 


Se = ^(tr(iy[yo(+oo)]) -tr(iy[yo(-oo)])) 


-4 


3(27rr)2i?6 


(t<%oo)%oo) - ^i-oo)^(-oo) 


3(27rr)2i?6 


(2.16) 




2=1 


In the following sections we will focus on the most elementary process in which a single 
membrane splits into two. This corresponds to an instanton solution in which is 
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taken to be the irreducible representation J and is the representation Ji © J 2 , where 

Ji + J 2 = J SiS required by angular momentum conservation. For this process the Euclidean 
action is 


Se = 


1 




2 y.. j 1 JJ1J2 

{2TTT)Ry ^ ^ 8 N 


(2.17) 


where in the final equality we have used the relations T = l/[(27r)^lp] and (R/lp)^ = 2 ^'k‘^N. 

The (anti-)instantons described above are local minima of the Euclidean action of the 
pp-wave matrix model. Their contribution to physical transition amplitudes associated with 
membrane splitting/joining processes can be evaluated using a standard semi-classical ap¬ 
proximation. The latter involves the integration over the bosonic and fermionic collective 
coordinates associated with zero modes in the instanton background. Erom the semi-classical 
calculation of transition amplitudes it should be possible to extract effective interaction ver¬ 
tices for membranes. In the present paper we focus on the construction of classical instanton 
solutions, leaving the semi-classical evaluation of the corresponding amplitudes for a future 
publication. Among the questions that we will not address is whether the BPS instantons 
discussed here give the dominant contributions to the amplitudes, or instead other saddle 
points associated with non-BPS configurations need to be taken into account. 


2.2 Some properties of BPS instantons and their moduli spaces 


A key observation used in the analysis of 18 is that (2.12) can be mapped to the Nahm 


equations arising in the construction of BPS multi-monopole configurations [19|. In order to 


reduce (2.12) to the Nahm equations one considers the change of variables [1^ 

Z* = C'e2*/^y^ 


s = e-2*/^, 


(2.18) 


where C is a constant, which we take to be real, so that the Z*’s are Hermitian matrices. 


Substituting into (2.12) gives 


d Z* 


(2.19) 


which by a suitable choice of the constant C can be brought to the canonical form of the 
Nahm equations. 

In terms of the original variables, we are interested in solutions T*(f) approaching constant 
configurations proportional to different SU(2) representations at t = ± 00 . After the change 


of variables (2.18) we therefore need to consider Nahm equations defined on a semi-infinite 


interval with boundary conditions 18 


Z\s) 


L) 


{-00) 


+ ••• for s —)• 00 , Z''{s) 


L) 


(+00) 


+ 


for s —)• 0 , 


( 2 . 20 ) 


where the ellipses indicate sub-leading terms. 

The Nahm equations have been extensively studied for their central role in the physics 
of monopoles and a great deal is known about their properties. The particular boundary 


conditions (2.20) relevant in the present case are non-standard and less studied, although 


they have been considered in a different context in 124 . In general, little is known about 
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explicit solutions to (2.19)-(2.20). The only known examples appear to be solutions presented 
in [18] . These include a particular conhguration interpolating between a vacuum associated 
with the representation trivial vacuum corresponding to the J-dimensional 

representation 1 © 1 © • • • © 1. The explicit form of this solution, in terms of the original 
matrices, is 


Y\t) = 


1 


L] 


{2txT)B? 1 + e2h-*o)/R (-“) 


( 2 . 21 ) 


In the following we will refer to (2.21) as the Bachas-Hoppe-Pioline (BHP) solution. The only 


other known explicit solutions are slight variations of this one, obtained by means of tensor 
products. They interpolate between different specihc pairs of representations and were also 
presented in [18] . However, no systematic approach to the construction of solutions has been 
proposed and no concrete examples are known for the class of processes we are interested 
in, i.e. those in which all the membranes involved in the splitting or joining transition carry 
large angular momenta. In the following sections we propose a strategy for the construction 
of instanton configurations in this large angular momentum sector. For the important case of 
the splitting of one membrane into two we also present an explicit solution, which represents 
an essential step in the calculation of the associated physical transition amplitude. 

Various properties of the moduli space of BPS instantons were discussed in , where, 
in particular, the dimension of the moduli space of a general instanton configuration was 
determined. For an instanton describing the splitting process with m membranes at t = —oo 
and n membranes at t = +oo - i.e. connecting the two representations Ji © J 2 © • • • © Jm 
and © J 2 © • • • © (with — Y17=i ~ ^ n > m and assuming that this is an 

allowed process according to the criterion in [^) - the (complex) dimension of the moduli 
space, Jt4, was found to be 


dime {M) = ^(2i - 1) J' - ^(2i - 1) J*, 


( 2 . 22 ) 


i=l 


i=l 


where in the sums the representations are ordered by decreasing dimension, i.e. Ji > J 2 > 
• • • > Jm and J'i> J' 2 > ■ ■ ■ > J'n- For the splitting of one membrane into two (2.22) gives a 
(real) dimension 

dim]R(AI) = 2(J(+3J2 -Ji) = 4J2, (2.23) 

where we used Ji = J[ + J 2 = J. 

An interesting feature of the moduli space of BPS instantons is an additivity rule |18j . 
Given three vacua associated with the SU(2) representations and T(c'), the moduli 

spaces M.[a,b)^ and M.(a,c) of configurations connecting A and B, B and C and A 

and C respectively satisfy 


dim {Mi^a,b)) + dim (M(b,c)) = dim {M(^a,c)) ■ 


(2.24) 


In section 4.3 we will comment on how this property is reflected in the description of the 


instanton moduli space arising from our reformulation in terms of the Laplace equation. 

The instanton conhgurations discussed in the previous subsection interpolate between 
vacua of the pp-wave matrix model containing different numbers of membranes. In it 
was shown, based on properties of the representations of the relevant supersymmetry algebra, 
that these vacua are non-perturbatively protected and their energies are exactly zero in the full 









quantum theory. This means that different vacua do not mix and tunnelling transitions should 
be possible only if excited states are involved. It is natural to expect that the mechanism 
forbidding transitions between vacua should be associated with selection rules induced by the 
integration over fermion zero modes in the instanton background. This was shown explicitly 


m 


17 in the case of the special solution (2.21) and similar selection rules should exist for 


more general tunnelling amplitudes such as those that we study in this paper. 

3 Mapping to three-dimensional Laplace equation 

In this section we present our general strategy for the construction of instanton solutions 
describing tunnelling configurations interpolating between vacua containing large membranes. 


3.1 Continuum approximation 

The central ingredient in the approach that we develop is a continuum approximation valid 
for large J, i.e. for large matrix size. In this approximation scheme, the matrices, are 

replaced by functions, u^), of two spatial coordinates, (u^, d^), and (Euclidean) time, 

t. The mathematics underlying this replacement is essentially the same that is utilised in the 
derivation of the matrix model as a regularisation of the supermembrane theory [1,26. In 
the matrix regularisation, the starting point is the continuum theory which is formulated in 
terms of functions on the membrane 2+1 dimensional world-volume. At fixed time, for any 
function, /, defined on the membrane world-space parametrised by the coordinates ((t^,(T^), 
one introduces a corresponding matrix, /?(/). The map p between functions and matrices 
satisfies the properties 


Pifd) ~ 2 


+ 


Pi{f,9}) 
r f /dV 


2ttJ 


1 


J 


[p{f),p{9)], 


(3.1) 

(3.2) 

(3.3) 


where the two sides are understood to be equal up to terms of higher order in 1/J. The Lie 
bracket, {/, 5 }, between two functions, /((T^,(T^) and is defined as 


{f,9} = 


df dg df dg 


da^ da‘^ da^ 

and the conventional constant [<t] denotes the area of the base space 


(3.4) 


kl = / dV. 


(3.5) 


In the present case we use a substitution of this type in reverse order. The starting point 
are time-dependent matrices, which we approximate by functions depending on two auxiliary 
continuous coordinates. Such an approach, based on the approximation of discrete objects by 


^For a discussion of the matrix regularisation emphasising its interpretation as an approximation between 


discrete and continuum objects, see 27 . 
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functions of continuous variables is of course quite standard. An elementary example of this 
type of approximation is the description of sound waves (phonons) in a solid, i.e. a discrete 
lattice, by a continuum wave equation. 


Using the rules (3.1)-(3.3), the instanton equations (2.12) obeyed by the matrices Y^{t) 
get mapped to non-linear partial differential equations for a set of functions 


, cr 


dy" 

dt 




{2ttT)R [cr] j r ,■ n 

dvr J 


(3.6) 


The boundary conditions for the matrices T*(t) at t = ±oo determine the corresponding 
boundary conditions for the functions ?/*(t, u^). The latter should describe families of 


concentric spheres at t = Too. The approximate identities (3.1)-(3.3) imply that, for large J, 


solutions to the original equations (2.12) are well approximated by functions satisfying these 


equations with the appropriate boundary conditions. 


The linear term in the continuum version (3.6) of the BPS equations can be eliminated 
by the change of variables 

.. (2vrr)i?2 


dvr 




s = e-2V«, 


(3.7) 


which is analogous to (2.18) that was used in 
equations 


Substituting into (3.6) yields the 


M ,wz,r 1 


(3.8) 


for the functions (s, , a"^), i = 1,2,3, where the variable s is dehned on a semi-infinite 
interval with s = -|-oo corresponding to t = —oo and s = 0 corresponding to t = -|-oo. As in 
the case of the y*’s the boundary conditions for the z*’s at s = -|-cx) and s = 0 follow from 


those required for the corresponding matrices (2.20). We will return to the analysis of the 


relevant boundary conditions for the equations (3.8) in the next section where we discuss 
various explicit solutions. 


The equations (3.8), sometimes referred to as the SU(oo) Nahm equations (after an ap¬ 


propriate rescaling to normalise the coefficients), are known to be (locally) equivalent to the 
three-dimensional Laplace equation 20 21 ^ In order to map (3.8) to the Laplace equation 
one uses a so-called hodograph transformation - a change of variables in which the roles of 
dependent and independent variables are exchanged, 




) ^(s,^7^C^2) 


(3.9) 


A fundamental feature of such a transformation is that it allows one to map a non-linear 
differential equation to a linear one. 


The change of variables (3.9) means, in particular, that one considers the originally in¬ 
dependent variable s as a function, i?i( 2 *), of the new independent variables z®, i = 1,2,3. 
One can then show that if the z®’s satisfy the equations (3.8) as functions of then 

s = (piz"^) obeys the Laplace equation 




(3.10) 


^ This change of variables is similar to that used in to relate the Poincare coordinates of the AdS space 
to the coordinates of the pp-wave space. It is also reminiscent of the coordinate transformation involved in 
the radial quantisation of the CFT on the boundary. 


The continuum Nahm equations were also considered in connection with membrane theory in 28 . 
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A geometric derivation of this result is presented in appendix [A) Alternative proofs can be 
found in (^[^ . 

This reformulation of the continuum version of the Nahm equations has a simple and 
very intuitive interpretation. The time variable s becomes a ‘potential’ function obeying the 
Laplace equation. The associated equipotential surfaces, i;i)( 2 ;*) = const., for different values of 
the potential represent implicit equations describing the membrane profile at the given time. 
Therefore a sequence of equipotential surfaces for values of s ranging from s = +00 to s = 0 
directly captures the evolution of the membrane configurations in Euclidean time. 


3.2 Classical membrane perspective 


In the previous subsection we have obtained the equations (3.6) as an approximation, valid for 


large matrices, to the original instanton equations (2.12). More specifically the approximation 


requires the dimension of all the irreducible representations contained in the vacua at t = ±00 
to be large. The same equations can also be understood as classical (Euclidean) equations of 
motion for the membrane theory in the AdS 4 x background in the pp-wave approximation. 
In the notations of the Euclidean action for the pp-wave membrane theory (in the Aq = 0 
gauge) is 


Se = / dtdV 


1 J 


dy^ 


1 J 


dx^ 


+ 


1 J 


2 J 




+ ^ 1)1 {{l”,l")" + 2 {x”,!,‘}" + /) 

3z J 1 rtT„,123/ 


1 J 




2 [d 


(3.11) 


As in the case of the matrix model, the bosonic part of this action can be rewritten as a sum 
of squares plus a boundary term, 

. 2 


[ dtd^cr 
2R d ./ 


dx"^ 


1 , ^.2 [d^ 

+ 




(3.12) 




2 

R' 


RT [d 


^ (1 

dt yR' 


RT id 


This formula shows that (3.6) can be obtained minimising the membrane Euclidean action 


and therefore these equations describe the BPS instanton configurations of the membrane 
theory. 

This is of course consistent and not surprising. The matrix model contains membrane 
degrees of freedom and the matrix configurations we focussed on represent regularised mem¬ 
brane states. However, we prefer to emphasise the point of view presented in the previous 


subsection in which the continuum equations (3.6) are viewed as an approximation to the 


corresponding matrix model instanton equations. This is because the matrix model itself is 
more fundamental as a candidate for a microscopic formulation of quantum M-theory. More¬ 
over the calculation of physical transition amplitudes in semi-classical approximation should 
be carried out in the matrix model. This in particular should allow one to compute tunnelling 
amplitudes between states corresponding to non-Abelian degrees of freedom associated with 
excitations in off-diagonal blocks, which have no simple counterpart in the continuum. 
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4 Solutions 


In this section we discuss solutions to the Laplace equation which correspond to membrane 
splitting processes. We will see that, in order to construct solutions with the required prop¬ 
erties, we need to introduce the concept of Riemann space. After explaining the general 


features of our proposal, we present exact solutions to the Laplace equation (3.10) with 


appropriate boundary conditions, which, based on the arguments in the previous section. 


provide approximate solutions to the BPS instanton equations (2.12). 


We first discuss the simplest solution corresponding to a static single spherical membrane, 
i.e. the simplest stable vacuum of the pp-wave matrix model. This case (which can be thought 
of as a solution in the zero instanton sector) allows us to illustrate the effects of the sequence 
of changes of variables that we use to map the original instanton equations to the Laplace 


equation. In section 4.2 we then discuss the solution corresponding to the splitting of a single 
membrane into two membranes, which requires the introduction of the notion of Riemann 
space. More general splitting processes, which involve more complex examples of Riemann 
spaces, are discussed in section 4.3 A reformulation of the solution (2.21) discussed in m 18 


in terms of the Laplace equation is presented, together with other simple new solutions, in 
appendix 

4.1 Stable sphere 

In order to understand the asymptotic behaviour of non-trivial instanton solutions and the 
corresponding required boundary conditions, it is instructive to first consider the simplest 


solution to the BPS instanton equations (2.12), namely the static configuration corresponding 
to the irreducible representation J, 


Y\t) = 


jL\ Vt G (—oo, -l-oo), 


(27rr)R2^ 

where [L^,lJ] = The continuum counterpart of this solution is 

y\t,cr^,a'^) =rn\a^,a‘^) , 


(4.1) 


(4.2) 


where r = J/[(2TrT)R‘^] and is a unit vector in the radial direction. Using the 

standard parametrisation of the sphere we can take 


n*(cr^, cr^) = (sin6lcos(/J, sin0sin</5, COS0), with cr^ = l — cos0, = (p . 


(4.3) 


Equations (^4.2l-(4.3) provide a static solution to the continuum version of the BPS instanton 


equations (^3.6^, as can be verified using the definition of the Lie bracket (3.4). In (4.3) we 


have chosen a parameterisation suitable for the North pole patch of the unit sphere. In 
general there is freedom in the choice of the (ci^, cr^) variables associated with the invariance 
of the membrane theory under area preserving diffeomorphisms. The (cr^,cj^) coordinates 
should be such that the area element dcr^du^ coincides with the natural SO(3) symmetric 
area element on a round sphere. 


Using the change of variables (3.7) we get a solution to the continuum Nahm equa- 

(4.4) 


tions (3.8) 


J 1 


z*(s,ct\cj^) = -- n \ a ^, a ^) 

47r s 
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This represents a sphere with radius increasing as s varies from +oo to 0. Inverting (4.4) to 


obtain s = we immediately deduce that the potential solving the associated Laplace 

equation satisfies 


47r|z| 


where z = 


Therefore the solution to the Laplace equation corresponding to a 
single static membrane conhguration is simply the Coulomb potential generated by a positive 
point charge J located at the origin. 

We note that a static spherical membrane in the (t, y) variables becomes a time-dependent 
spherical membrane in the (s,z) variables, with radius changing in time as 1/s. Using the 


transformation (3.7) in reverse order, one finds indeed that (4.5) corresponds to a membrane 


whose distance from the origin is constant in time t, 


J 


\y\ = 


{2'kT)E? ’ 


with = 


(4.6) 


which of course represents a static sphere of radius J/{2ttTB?). 


As a simple generalisation of (4.5) we can consider the case of a point charge located away 
from the origin. The potential is 


(/{z) = 


J 


47r|z — ZqI 


Zo / 0. 


(4.7) 


In this case the equipotential surfaces are spheres centred at zq. 
variables we have 

I I _ J 1^1 

“ (27rr)i?2 Iz-zol' 


In terms of the original y 


(4.8) 


This configuration corresponds to a spherical membrane moving in towards the origin from 
infinity. As (Euclidean) time evolves from t = —oo, corresponding to s = +oo where |z — 
ZqI —)• 0, to t = + 00 , corresponding to s = 0, \y\ decreases from 00 to a constant. This 
run-away behaviour of the solution in the infinite past makes it non-physical, because the 
solution does not satisfy the required boundary conditions. This is a general result, which 
implies that the potentials relevant for the description of configurations interpolating between 
stable vacua of the pp-wave matrix model can only involve positive charges located at the 
origin 


4.2 Single membrane splitting and Riemann space 

In this section we consider the most interesting solution, which describes a single membrane 
with angular momentum J splitting into two membranes with angular momenta Ji and J 2 
(with Ji -|- J 2 = J). This corresponds to a configuration interpolating between a spherical 
membrane of radius J/{2ttTR'^) and two concentric membranes of radii Ji/{2 'kTR?), i = 1, 2. 
It is the most elementary example of a splitting process and it allows us to illustrate general 
properties which are common to all splitting/joining transitions. 

^ The location of negative charges can be arbitrary. In appendix we will present explicit examples of 
physically acceptable solutions involving negative charges away from the origin. 
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4.2.1 Qualitative picture: introduction of the concept of Riemann space 

We will now argue that it is not possible to describe a membrane splitting process in terms 
of a solution to the Laplace equation formulated in a standard Euclidean space Instead, 
we are led to introduce the concept of Riemann space explained in detail below. In order 
to motivate our proposal, let us first discuss the qualitative features of a typical instanton 
solution describing a splitting process. The transition from a spherical membrane to two 
concentric membranes can be qualitatively described as involving the steps depicted in fig¬ 
ure]^ Starting with a single membrane of radius J/{2ttTB?) a small perturbation forms and 
grows until the membrane splits into two, resulting in the creation of a second disconnected 
membrane inside the first one, as shown in steps (i) to {in) in the figure]^ 

To understand the subsequent steps it is important to notice that membranes are charged 
and thus oriented objects. In figure the orientation of the membranes is indicated by 
the relative position of the continuous and dashed contour lines: membranes have an outer 
continuous contour, whereas in anti-membranes the outer contour is dashed. In the pp-wave 
background a configuration of two concentric membranes with the same orientation, which 
is BPS, is stable as a result of the balance between the gravitational attraction and the 
force associated with the three-form flux. The two forces add up in the case of a membrane 
anti-membrane pair, resulting in an unstable configuration. 

As shown in figure after the splitting takes place the internal membrane that is formed 
has opposite orientation compared to the initial one and the resulting configuration is not 
stable. To complete the transition to a stable vacuum consisting of two concentric mem¬ 
branes it is necessary for the internal membrane to flip its orientation. This process requires 
intermediate steps in which the membrane self-intersects after developing cusps, see {iv)-{vi) 
in figure 

The process shown in the figure has axial symmetry, but this is of course not a require¬ 
ment. In general the details of the intermediate configurations can vary, however, an essential 
element, common to all splitting processes, is the fact that a flipping of the orientation is 
inevitable. In the illustrative example shown in the figure the sequence of steps leading to the 
final two-membrane state involves a splitting resulting in the formation of a second internal 
anti-membrane followed by the flipping of the orientation of the latter. However, as we will 
see in explicit examples in the next subsection, the splitting and flipping steps can also take 
place in reverse order, or simultaneously. 

Applying the approach described in section]^ our goal is to construct a solution, 4>{z), to 
the Laplace equation such that the associated equipotential surfaces reproduce the qualitative 
behaviour shown in figure 

The crucial feature of any splitting process, which emerges from the qualitative arguments 
we presented above, is the fact that necessarily there exists a finite interval of values of the 
potential for which the equipotential surfaces are self-intersecting. Such a behaviour is not 
compatible with the potential (i){z) being a solution to the Laplace equation in the ordinary 
]R3 

space, because the associated gradient, V(/)(z), would be ill-defined, since the normal 
direction differs on the two portions of a self-intersecting equipotential surface. In the case of 

® The point where the splitting occurs is singular from the point of view of the partial differential equation 
governing the evolution of the membranes. The nature of the singularity and the properties of the instanton 
solution in the vicinity of a generic splitting point are universal and can be studied using the description in 
terms of the Laplace equation, as discussed in appendix |B| 
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(v) 


(vi) 


Figure 1: Qualitative depiction of the splitting process. The process necessarily involves both 
splitting and flipping. 


the intersections occurring during the flipping of the membrane orientation, this issue arises 
not just at isolated singular points, but over a finite region corresponding to an interval of 
values of (f){z). This motivates us to consider multi-valued potential functions. 

More concretely the above observations lead us to propose that the potential relevant 
for the representation of a membrane splitting process should be a solution to the Laplace 
equation in a space, which we refer to as a Riemann space, that is a three-dimensional 
generalisation of a two-sheeted Riemann surface. Such a space consists of two copies of 
connected by a branch surface, bounded by a branch curve or loop. For simplicity we will 
assume the surface connecting the two IR^’s to always have the topology of a disk and in 
the following we will use the expression branch disk without necessarily implying a circular 
planar shape Deformations of the branch disk which leave the branch loop fixed yield 
physically equivalent Riemann spaces. This is analogous to properties of branch cuts and 
branch points of Riemann surfaces. A space defined in this way is locally equivalent to IR^, 
but differs globally. The use of a Riemann space makes it possible to have a potential which 
locally solves the Laplace equation, while avoiding the problems caused by self-intersecting 
equipotential surfaces. The solution is such that different portions of self-intersecting surfaces 

® In appendix we present a reformulation of the Laplace equation in a Riemann space as an integral 
equation and we discuss in more detail the boundary conditions at the branch disk. To construct solutions 
relevant for the description of membrane splitting processes we require the potential to be finite at the branch 
loop and also to decay sufficiently rapidly at infinity. 


15 

























live on different sheets of the Riemann space, so that no intersections take place within the 
same copy of IR^. 


The description of a stable spherical membrane in section 4.1 provides an indication of 
what the asymptotic behaviour of the potential should be in the case of a splitting process. 
The infinite past, t —>■ —oo, corresponds to s —)■ +oo, i.e. large values of the potential. In this 
region the equipotential surfaces should be small spheres with radius proportional to J and 
approaching zero. The inhnite future, t —>■ +oo, corresponds to s —>■ 0, i.e. small values of 
the potential approaching zero from above. In this region the equipotential surfaces should 
be two large concentric spheres with radii proportional to Ji and J 2 and growing indefinitely. 
Using a Riemann space we can construct a potential which provides a concrete realisation of 
this type of asymptotic behaviour and is also well-defined in the intermediate region where 
the membranes self-intersect. More precisely, in order to account for the behaviour at large (j) 
we consider a point charge J at the origin of the first IR^, so that close to the charge we have 
a Coulomb potential, whose equipotential surfaces are small spheres. The behaviour in the 
two IR^’s for small s, on the other hand, can be understood as resulting from the splitting of 
the initial flux J of V(/>, with a fraction Ji going to infinity in the hrst space and a fraction 
J 2 = J — Ji passing into the second space through the branch disk. 

Figure [^illustrates qualitatively the splitting process as rendered in the Riemann space, 
with the six rows in the picture corresponding to the same stages of the splitting process shown 
in figure The first two columns depict the two sheets of the Riemann space, with each copy 
of IR^ represented as a rectangle. The horizontal slit corresponds to the branch disk, with the 
dots at the end points indicating the branch loop. The third column shows the two copies of 
IR^ superposed. Membranes, or portions of membranes, are represented as continuous lines 
in the first IR^ and as dashed lines in the second IR^. The same notation will be used in the 
figures throughout the paper. Re-analysing the splitting process as presented in the figure 
allows us to clarify how the issues associated with the flipping of the membrane orientation 
and the self-intersections in the equipotential surfaces can be addressed using a Riemann 
space. As the value of the potential decreases, away from the point charge, the equipotential 
surfaces intersect the branch disk and therefore extend partially into the second copy of IR^. 
This is shown at step {in), where the splitting has already taken place. Crucially the use of a 
Riemann space allows us to have self-intersecting surfaces which result from the superposition 
of portions of surfaces without self-intersections in each sheet. This is illustrated in row {v) 
in figure Notice that in the third column intersections occur only between a continuous 
and a dashed line, never between two continuous or two dashed lines. The final conhguration 
in the splitting process is represented by two membranes, each contained entirely in one copy 
of IR^, with the superposition resulting in two concentric surfaces. 

As already pointed out, a feature of the solution is that during the flipping process the 
equipotential surfaces develop cusps. These special points occur as the equipotential surfaces 
intersect the branch loop, as shown at step (iv) in hgure [^ From the explicit solutions 
discussed in the next subsection one can verify that the potential itself is regular over the 
entire Riemann space (except for the origin of the first IR^ where the point charge is located), 
including the branch loop where the cusps arise [^ At these points, however, the gradient 
of the potential diverges. This allows us to physically characterise the branch loop in the 
Riemann space in terms of the behaviour of the solution to the continuum Nahm equations. 


'This is true in general and not only in the cases with axial symmetry for which we present plots. 
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Figure 2: Qualitative depiction of splitting in Riemann space. 
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Recalling that the potential, s = (p{z), is obtained inverting the functions z'^{s), we deduce 
that the branch loop, where ^4>{z) diverges, corresponds to points where the velocity of the 
front of the membrane, as described by z^{s) in Euclidean time, vanishes. 

Various other properties of the instanton solutions corresponding to membrane splitting 
processes and their moduli spaces can be given an intuitive interpretation using a description 
in terms of Riemann spaces. We will discuss some of these aspects in section 4.3.2[ 


4.2.2 Analytic solution 

In section we have shown that the problem of solving the continuous version of the instan¬ 
ton equations can be mapped to that of finding solutions to the three-dimensional Laplace 
equation. Then in the previous subsection we have discussed the boundary conditions that 
we propose to consider in order to obtain a potential with the required properties. We now 
present analytic examples showing explicitly that solutions obeying such boundary conditions 
can be constructed. 


Based on the examples of static solutions presented in section 4.1 we expect the potential 
for large s (corresponding to t —>• — oo in the original Euclidean time variable) to have a 
Coulomb-like singularity. The considerations in the previous subsection about the expected 
qualitative behaviour of the potential then lead us to consider the Laplace equation in a 
Riemann space made of two copies of IR^, with boundary conditions associated with the 
presence of a single positive point charge at the origin of the first 

The idea of studying the Laplace equation in a Riemann space has actually been con¬ 
sidered long ago by Sommerfeld in a 1896 paper 23 , which introduces the idea of Riemann 


spaces to develop a generalisation of the standard method of images for the solution of elec¬ 
trostatics problems. The examples studied in 23 involve multiple copies of with branch 


curves consisting of straight lines. Eollowing Sommerfeld’s proposal, an analytic solution to 
the Laplace equation with boundary conditions precisely of the type we are interested in 
was constructed by Hobson in . This paper considers a Riemann space consisting of two 
copies of R^ connected by a flat circular branch disk and computes the potential generated 
by a point charge located in the first R^. More recently, Riemann spaces have been consid¬ 
ered in [30,31 . It can be shown that Sommerfeld’s solution with a straight branch line and 


Hobson’s solution with a circular branch loop can be related by an inversion transformation, 
see appendix [D| 

Hobson’s solution is written in terms of so-called peripolar coordinates, which are par¬ 
ticularly suited to the specific geometry under consideration. The peripolar coordinates of a 
point P are designated by {p, 9, ip). Their definition in R^ and their relation to the Cartesian 
coordinates, {z^,z^,z^), can be given as follows. One starts with a circle of radius a, which 
will later be identified with the branch loop and we conventionally take to lie in the 
plane. A plane containing the point P and the axis of the circle, which we identify with the 
z^ axis, intersects the circle at two diametrically opposite points, A and B. Denoting the 
distances of P from A and B respectively by r and r', the coordinate p of P is defined as 


= log - , 
r 


(4.9) 


while 9 is the APB angle. The angle ip is the standard polar angle for the projection, N, of 
P onto the {z^,z‘^) plane. Figure [^illustrates the definition of p, 9 and ip and their relation 
to the Cartesian coordinates. 
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Figure 3: Relation between peripolar and Cartesian coordinates 


The angle (p is taken to vary between 0 and 2ti. From the definition p G (— 00 , + 00 ) and 
6 is defined to be in the interval [—vr,+7r]. The angle 6 goes to 0 when the distance of P 
from the origin goes to 00 and in the region of the z^) plane outside the circle. It has a 
discontinuity as one passes through the interior of the disk bounded by the circle of radius 
a. It approaches +7r or —vr if one approaches a point inside the disk from above > 0) 
or below (z^ < 0) respectively. Points on the z^ axis correspond to p = 0 and points on the 
circle of radius a in the {z^, z‘^) plane have p = ± 00 . 

Denoting by ^ the distance of the point N from the origin, the peripolar and Cartesian 
coordinates are related by 


z^ = ^ cos , z^ = ^ sin ip , 


z3 = 


asin0 


cosh p — cos 0 


(4.10) 


where 




a sinh p 
cosh p — cos 9 


(4.11) 


To describe (in peripolar coordinates) a Riemann space consisting of two copies of con¬ 
nected by a branch disk coinciding with the disk bounded by the circle of radius a in the 
(z^, z^) plane, one simply allows the range of 9 to extend from —vr to Svr. The intervals [—vr, vr] 
and [vr, 3vr] correspond to the first and second IR^ respectively. Moving in the first IR^ from 
infinity towards the origin along the positive z^ axis (or any other curve in the 2 ;^ > 0 region) 
9 is positive and increases from 0 to vr. Crossing the disk at 2 ;^ = 0, one crosses into the 
second IR^, while 9 varies continuously. Moving along the negative z^ away from the origin in 
the second space, 9 continues to increase from vr, reaching 2vr as z^ —>• — 00 . Approaching the 
origin along the positive z^ axis in the second space the angle 9 varies from 27r to 3vr, which 
is reached on the upper side of the disk. Passing through the disk one crosses back into the 
first space (with z^ < 0) and 9 goes back to —vr. 
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The Laplacian in peripolar coordinates takes the form 


/ cosh p — cos 6 

V ^2 


, , /IX (1 — cos 0 cosh p) d 

(cosh p - cos 6) —^ H-—-TT 

ap^ sinh p op 


, , /i\ ^ (cosh p — cos 6) 92 

+{cosh - cos 0)^ - sm 9- + ^ 


(4.12) 


In Hobson computed the potential solving the Laplace equation in the two-sheeted Rie- 
mann space with a circular branch loop for an arbitrary relative position of the point charge 
in the hrst space relative to the branch disk. Denoting by zq the location of the point 
charge, with peripolar coordinates (po,0o,(po)) the potential at a generic point z of coordi¬ 
nates (p, 9, <p) is 


(/>(z,zo) = 


J 


47r|z — zqI 


11 . / ( 9-00 

- -|— arcsm cos —-— 

2 TT \ \ 2 


cosh a + 1 


where 

and 


cosh a = cosh p cosh po — cos((p — <po) sinh p sinh po 


|Z — ZqI = 


cosh a — cos(0 — 0o) 

(cosh p — cos 0) (cosh po — cos 9q) 


(4.13) 

(4.14) 

(4.15) 


Notice that the potential (4.13) is periodic in 0 with period dvr, as appropriate for the two- 


sheeted Riemann space. One can verify that (j){z,zo) is well-dehned and finite everywhere, 
including the branch loop, except for the location of the point charge, i.e. z = zq in the first 
space, where it has a Coulomb-like divergence. 

As anticipated in the previous subsection, the flux of V(() provides a measure of the 
angular momentum carried by the different membranes, which is also proportional to their 
respective radius. In the case under consideration we have a two-sheeted Riemann space 
with one point charge in the hrst copy of We will see that the equipotential surfaces for 
large values of the potential are small spheres centred at the location of the point charge. 
They represent the single membrane with angular momentum J at large negative values of 
the original Euclidean time t. The hux of through these surfaces is J. The asymptotic 
equipotential surfaces for small values of the potential approaching zero will be shown to be 
two (approximately) concentric spheres with diverging radii, one in each copy of IR^. The hux 
of V(/> through these spheres (coinciding with the hux at inhnity in the respective IR^) equals 
the angular momentum of the corresponding membrane, Jj, i = 1,2. Conservation of the 
hux, which corresponds to conservation of angular momentum for the membranes, implies 
J = Ji -|- J 2 . The way in which the hux J is split between the two spaces is controlled by 
the size of the branch disk and by its position relative to the point charge. For a disk of 
given radius, if the point charge is located very far from the disk only a small fraction of the 
hux passes through the disk into the second space and thus we have J 2 Ji, corresponding 
to a hnal state with one membrane much larger than the other. If we reduce the distance 
between the point charge and the branch disk, the fraction of hux passing into the second 
space increases. When the distance of the charge from the disk tends to zero, we approach 
the case in which the hux is equally spit between the two spaces, which in turn corresponds 
to having two membranes of equal radius in the hnal state. 
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These general considerations can be made more precise by studying the asymptotic be¬ 
haviour of (I){z,zq) at a large distance from both the point charge and the branch disk in 
each of the two spaces. For this purpose we analyse (l){z,zo) in the region defined by /? —)■ 0, 
0 —)• 0 in the first space and p —)• 0, 0 —27r in the second space. In both cases we have 
cosh a ~ coshpo- The asymptotic behaviour at infinity in the first space is 


(j){z,zo) ! 

and in the second space it is 
(j){z,Zo) s 


J 


47r|2: — Zq\ 


J 

47r|2: — ZqI 


1 1 

—I— arcsin 
2 TT 


1 1 

-arcsin 

2 TT 


COS ( ^ 


cosh (^) 


cos ( ^ 


cosh (^) 


(4.16) 


(4.17) 


The factors in square brackets in (4.16) and (4.17), multiplying the Coulomb potential 


J/{Att\z — zo|), represent the fractions of the total flux staying in the first space or escaping 
into the second space, respectively. They determine the angular momenta Ji and J 2 of the 
two membranes in the final state as fractions of the total angular momentum J. 

In the axially symmetric case in which the point charge is located on the axis of the 


branch disk, corresponding to po = 0, the asymptotic formulae (4.16) and (4.17) simplify and 
we get 


(j){z,zo) 


J 


47r|z — ZqI 


1 - 


N 

27r 


in the first space and 


(^(z,Zo) 


J 


l^ol 


(4.18) 


(4.19) 


47r|z — ZqI 27r ’ 

in the second space. Notice that for 9q —?• ±7r, i.e. when the location of the point charge 
approaches the disk, the ratio goes to 1/2 as it should be. 

In the following we set a = 1 and we present plots of the equipotential surfaces of the 


potential (4.13) in the axially symmetric case, po = 0. We can also set po = 0 without loss 
of generality, so that the only remaining parameter in the solution is the angle Oq, whose 
value controls the distance of the branch loop from the point charge As will be shown 
in the explicit examples below, the value of Oq also affects the order in which the splitting 
and flipping of the membrane orientation take place. Therefore the sequence in which these 
steps occur is correlated with the way the angular momentum gets divided between the two 
membranes in the final state. 

Using the relations ( 4.10| )-(4.11 ) in (4.13), one can obtain the form of the potential in 
Cartesian coordinates, cj) = (/>(z^, z^), which is the expression used to produce the plots 

presented below. In all the following figures we use a rescaled potential, (p, related to (4.13) 
hy (f)= (4:71/J)(j) 


In figure]^ we show contour plots for the potential (4.13) with Oq = —47r/5. The figure 
shows families of equipotential surfaces separately in the two copies of Each is 
represented by a square, with the vertical direction being the direction of the z^ axis. The 


“We note that solutions with different do can be related by a conformal transformation. 
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branch disk is indicated by a horizontal slit, the point charge (contained in the first which 
is on the left) is denoted by a small crossed circle below the disk. The contours are displayed 
as continuous lines in the first space and as dashed lines in the second space. 



Figure 4: Contour plots showing families of equipotential surfaces in each of the two copies 
of which constitute the Riemann space. The contours shown are for 6q = —^ and 
correspond to unequally spaced values of (j). The values of the potential used here are 
^ = {3.113,1.613,0.878,0.593,0.368,0.323,0.310,0.278,0.173,0.128} and include those used 
in figure 


The following figures depict the evolution of the profile of the membranes throughout 
the splitting process for different values of Oq. They show how the equipotential surfaces 
for the potential in the analytic solution (4.13) reproduce the steps that were qualitatively 
discussed in the previous subsection. In these figures the two copies of R^ are superposed. 
The (portions of) membranes living in the first or second space are depicted as continuous 
or dashed lines respectively. The figures display axially symmetric solutions, therefore the 
three-dimensional shape of the membranes can be generated rotating the contours about the 
vertical axis. 

It is interesting to notice that the equipotential surfaces deviate significantly from a 
spherical shape only for a rather small range of the potential around the value where the 
splitting takes place. This is the region in which the equipotential surfaces cross the branch 
disk. The surfaces then quickly revert to Coulomb-like behaviour outside this region for both 
larger and smaller values of (j). 

Figure]^ shows the membrane profiles corresponding to the choice Oq = —Air/b in (4.13). 
In this case the point charge is quite close to the branch disk and correspondingly the 
flux/angular momentum gets split almost evenly between the two spaces. The sequence 
of plots shows the formation of cusps, which in these examples with axial symmetry occur 
simultaneously at all points on the branch loop. The membrane is then seen self-intersecting 
before it splits into two. As a result, with this choice of boundary conditions when the split¬ 
ting occurs the second membrane has already the correct orientation, but the two membranes 
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are still intersecting. Notice that, as anticipated in the qualitative description of the previous 
subsection, all the intersections between equipotential surfaces in the figure (in the third, 
fourth and fifth panel) involve a dashed and a continuous line, i.e. they always occur between 
(portions of) equipotential surfaces belonging to different sheets of the Riemann space. The 
same feature can be observed in all the subsequent figures in this and the next subsection. 

Figure]^ shows the evolution of the membrane profiles for Oq = — tt/G, which corresponds 
to a point charge much further from the branch disk. In this case the splitting process 
follows a sequence of steps similar to those depicted qualitatively in figures and The 
initial membrane splits into a membrane and an anti-membrane, which subsequently flips 
its orientation. The cusps are formed after the splitting and the subsequent flipping of the 
orientation of the internal membrane requires the latter to self-intersect. This takes place in 
the fifth panel in the figure and is shown more clearly in figure which shows an enlargement 
of the rectangular area marked in figure This example shows a concrete realisation, in an 
exact analytic solution, of the crucial step - illustrated at stage (v) in figures and [^- which 
led us to the concept of Riemann space. Notice the relative size of the two membranes in the 
final state; in this case (with the point charge further away from the disk) the flux is split less 
evenly than in the previous case and correspondingly the second membrane is much smaller. 

Figure shows the splitting process for 0o = —vr/3, i.e. a value intermediate between 
the cases shown in the previous two figures. In this case the splitting and flipping phases 
occur simultaneously. At the splitting point the second membrane is still self-intersecting and 
its orientation not completely flipped. The relative size of the two membranes in the final 
state reflects the fact that in this case the value of Og (which controls the distance between 
point charge and branch disk) is intermediate between those in the two previous examples. 
It corresponds to a splitting of the flux/angular momentum J into Ji and J 2 such that the 
difference | Ji — J 2 I is larger than in the first case and smaller than in the second case. 

If the point charge is not located on the axis of the disk the qualitative behaviour of the 
solutions and the associated equipotential surfaces remains the same. The processes display 
similar sequences of steps involving splitting and flipping. In this case, however, the cusp 
singularities where V(/) diverges do not form simultaneously on the whole branch loop, but 
rather at pairs of points. 

The equipotential surfaces for the potential (4.13), once expressed in Cartesian coordi¬ 


nates, represent implicit equations of the membrane profiles as described by the z^{s) vari¬ 
ables. The evolution of the shape of the membranes during a splitting process in terms of 
the original coordinates, can be reconstructed using the change of variables (3.7) in the 


potential 4>{z^, z^) The qualitative features of the equipotential surfaces parameterised 

by the y* coordinates are similar. The only significant difference is that the size of the mem¬ 
branes does not grow steadily as seen in the previous figures. This is the same difference that 
was observed in section |4.1| in the case of the description of the static spherical solution in 
terms of the y* and 2 * variables. Using the y* variables one can instead verify immediately 
that the solutions describe a process in which the final membranes have radii which add up 
to the radius of the single membrane in the initial state. In the next subsection we will 
present plots showing the evolution of the membrane profiles, using the parametrisation in 
terms of the y* variables, in the case of solutions interpolating between configurations with 


®In addition one also needs to redefine the j/* variables by constant shifts to bring the location of point 


charge to the origin and avoid the run-away behaviour mentioned at the end of section 4.1 


23 








$=3,113 




$=0,278 



Figure 5: Evolution of the membrane profiles, showing the different phases of the splitting 
process in Euclidean time in the case Oq = 
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Figure 6: Evolution of the membrane profiles, showing the different phases of the splitting 
process in Euclidean time in the case Oq = —7r/6 
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Figure 7: Enlargement of the area marked by a rectangle in the fifth panel of figure 6 


two membranes in the initial state and two in the final state. 


4.3 General Riemann spaces 

In the previous subsections we presented explicit solutions to the Laplace equation describing 
the most elementary splitting process with a single membrane in the initial state and two 
membranes in the final state. It is natural to ask whether a similar approach can be used 
to construct other solutions corresponding to (anti-)instantons of the pp-wave matrix model 
interpolating between states with different numbers of membranes. A large class of solutions 
connecting multi-membrane states was shown to exist in [18] . For our reformulation in terms 
of the Laplace equation to be equivalent to the original (continuum) instanton equations, 
solutions should exist for all the allowed instanton configurations. Therefore it should be 
possible to identify suitable boundary conditions corresponding to all the pairs of initial and 
final states for the original equations, which are permitted according to the criteria in |18| . 
Indeed, we will show that the conditions for the existence of BPS instanton solutions obtained 


m 


18 can be reproduced using our formulation in terms of the Laplace equation in a Riemann 


space. 

An analysis of the general features of the solutions we discussed above provides hints for 


the generalisation to more complicated processes. In the cases discussed in section 4.2 the two 
membranes in the final state correspond to equipotential surfaces in the asymptotic regions 
at infinity, one in each of the two copies of IR^ constituting the Riemann space. On the other 
hand the single membrane in the infinite past corresponds to small spheres centred at the 
location of the point charge. These considerations lead us to propose the following general 
prescription. In order to describe a process with n membranes in the final state one should 
consider a Riemann space made of n copies of Moreover, if the initial state contains m 
membranes one should consider m (positive) point charges at the origin of m distinct IR^’s. 
The values of the point charges, Ji, i = 1,2,... ,m, correspond to the angular momenta of the 
m membranes in the initial state. The values of the flux of V(^ at infinity in each copy of IR^, 
J', i = 1,2,... ,n, correspond to the angular momenta of the n membranes in the final state. 
We point out that an immediate consequence of this prescription is that by construction the 
number, n, of membranes in the final state is greater than or equal to the number, m, of 
membranes in the initial state. Therefore the necessary condition, n > m, for the existence 
of instanton solutions is automatically satisfied. 

In the case of a generic Riemann space the relation between the outgoing fluxes at infinity. 
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Figure 8: Evolution of the membrane profiles, showing the different phases of the splitting 
process in Euclidean time in the case Oq = —tt/S 
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J', and the values of the point charges at the origin, Jj, is non-trivial. This is because the 
way in which the flux of V(/> gets split among the different spaces depends in a complicated 
way on the number, shape and size of the branch disks and on their positions relative to the 


point charges. In the most elementary case of the solution (4.13) we gave explicit formulae for 


the asymptotic fluxes associated with the two membranes in the final state in (4.16)-(4.19). 


General multi-sheeted Riemann spaces can involve different combinations of branch sur¬ 
faces. This is the case even in the most elementary example discussed in the previous sub¬ 
section, in which one membrane with angular momentum J splits into two membranes with 
angular momenta Ji and J 2 . We described this process using a Riemann space made of two 
copies of IR^ with a single point charge, J, and we presented an exact solution involving a 
circular branch disk. However, the same process receives contributions from more compli¬ 
cated Riemann spaces in which the two IR^’s are connected by multiple branch disks, so long 
as the flux at infinity in the two spaces is split into the same fractions, Ji and J 2 . Of course 
in order to compute a specific physical transition amplitude it is necessary to sum all the 
contributions with the given initial and final states. The issue of how the different copies of 
1R3 

are connected is related to the parameterisation of the moduli spaces of the associated 


instanton solutions and we will briefly comment on this point in section 4.3.2 


At the present stage we do not have a proof of existence of solutions to the three- 
dimensional Laplace equation in arbitrary Riemann spaces, but we propose as a conjecture 
that solutions should exist for all the boundary conditions which according to our prescrip¬ 
tion correspond to allowed (anti-)instantons. The results of [23] and [29] provide a starting 
point for the analysis of the issue of existence of solutions in general Riemann spaces. One 
may view the explicit solutions in [23| and [29| as playing a role similar to that played by the 
explicit solution for a spherical conductor obtained by the method of images in connection 
with the theory of the electric potential for arbitrarily shaped conductors. Moreover the 
Laplace equation in the Riemann space can be recast into the form of an integral equation 
as discussed in appendix |E| This reformulation may also provide a useful approach to obtain 
a proof of existence for the solutions. 


4.3.1 More solutions 

Before discussing how some general properties of the instanton moduli space arise within 
our formulation using the Laplace equation, we now present other exact solutions. We use 
these specific examples to illustrate some interesting features arising more generally in the 
description of instanton configurations in terms of equipotential surfaces. 

The linearity of the Laplace equation of course implies that linear combinations of so¬ 
lutions are also solutions. We can take advantage of this simple observation to obtain new 
interesting examples. In particular, a solution describing a process with two membranes with 
angular momenta Ji and J 2 in the initial state and two in the final state - which based on 
our prescription requires a two-sheeted Riemann space with a point charge in each copy of 


IR'* - can be obtained as a linear combination of the potential (4.13) with charge Ji and 


the analogous potential with point charge J 2 at z = zq in the second space. To obtain the 
potential due to a point charge at the point corresponding to zq in the second space, we recall 
that in peripolar coordinates the first sheet is parameterised by 0 G [—vr, tt] and the second 
sheet by 0 G [tt, 37r]. Therefore the potential induced by a charge located at the point in the 
second IR^ corresponding to zq is obtained substituting 60 —)■ 0o + 27r in (4.13). The resulting 
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solution in the two sheeted Riemann space with point charges in both spaces is 


(l){z,Zo) = 


Jl 


+ 


47r|2: — Zo\ 

J 2 

ArK\z — Zq\ 
1 

47r|z — Zq\ 


11 . / fe-eo' 

—I — arcsm cos - , ,, 

2 vr \ \ 2 / V cosh a + 1 


1 1 . fe-9o 

-arcsm cos - 

2 vr I V 2 


{J1 + J2) , {J1-J2) 

-^-arcsm 

2 TT 



(4.20) 


cosh a + 1 


where \z — zqI ^md cosh a are given in (4.15) and (4.14) respectively and we have chosen the 


radius of the branch loop to be a = 1. The sign difference between the first two lines in (4.20) 
comes from the shift in 9q. Recalling that cosh a is independent of 9q, we get 


arcsm 


cos 


9 ~ { 9 q + 27r) 


cosh a + 1 


= — arcsm 


cos 


9-9o 


cosh a + 1 


(4.21) 

because the shift by vr in the argument flips the sign of the cosine and the inverse sine 


function is odd. Notice that if the two charges are equal, Ji = J 2 , (4.20) reduces to a simple 


Coulomb potential in both sheets of the Riemann space. This is consistent with our general 
prescription. With equal charges the fluxes of V(/> from the first space into the second and 
from the second space into the hrst are equal, irrespective of the position, zq, of the charges. 
So the flux at infinity in the two spaces is the same and hence the final state is the same 
as the initial state. Therefore this case corresponds simply to a stable vacuum conhguration 


with two membranes of the same size. Note that, up to an overall scale, the potential (4.20) 


can also be obtained from the superposition of the solution (4.13) and a simple Coulomb 
potential. 

We now focus on the special case of axially symmetric solutions (corresponding to po = 


Pq = 0) with Jl = 2J, J 2 = J for which (4.20) becomes 


4>{z,zo) = 


J 


47r|z — ZqI 


3 1 . / f9-9o 

- H— arcsm cos —-— 
2 TT y y 2 


cosh /? + 1 


(4.22) 


where we used the fact that cosh a. = cosh p when zq is on the axis of the disk. 


The final states in the process described by the potential (4.22) depend on the way in 


which the total flux at infinity is divided between the two spaces, which in turn is controlled 
by the value of 9q. Different choices for 9q give rise to solutions displaying various interesting 
features and below we present a few selected examples. 

The following figures HKH] and show the evolution of the two membranes described 


by the solution (4.22) for 9q = 9q = —tt/IO and 9q = —vr/G respectively. As in 


the previous subsection, the plots show equipotential surfaces for the rescaled potential (j) = 
(47r/J)i?i). The figures depict the equipotential surfaces for , y'^, y^) expressed in terms of 


the coordinates, y*(t), in the original instanton equations (3.6). The conventions used are the 
same as in previous figures. The (portions of) membranes in the hrst copy of IR^ are shown as 
continuous contours and those in the second as dashed contours. Since here we are using 
the y* coordinates, the contours represent the evolution of the prohles of the membranes from 
t = —00 to t = + 00 . Notice that with this parameterisation of the solution the membranes do 
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not expand steadily. However, their relative radii change between the initial and final state 
because of the transfer of angular momentum. The use of these variables makes manifest the 
conservation of angular momentum (JJ + J 2 = Ji + ^ 2 ); which in the figures is reflected in 
the fact that the radii of the spheres in the initial and final states add up to the same total. 
In the hgures we have not displayed the branch loop as in these coordinates its position is 
not constant. 

In the case shown in hgure the value of 9q = and thus the point charges are 

relatively close to the branch disk. As a result there is a large fraction of flux passing from 
the first into the second. This corresponds to a large transfer of angular momentum and 
thus a significant change in the relative radii of the two membranes. In the sequence shown 
in the figure the second and third steps show both membranes extending across both spaces. 
In the third panel the outer membrane develops a self-intersection. The small region marked 
by a square is shown enlarged in figure [TOl where the self-intersection is more clearly visible. 
The membranes then split and reconnect in the way shown in the fourth panel of figure 
resulting in two equipotential surfaces each contained entirely in one copy of IR^. 

In the case shown in figure [IT] the point charges are at 9q = —vr/lO, i.e. further from the 
branch disk and this means that there is a smaller transfer of angular momentum. With this 
choice of boundary conditions we observe another interesting feature: for some intermediate 
values of Euclidean time there are three membranes involved in the process. A small mem¬ 
brane detaches from the larger membrane in the first and subsequently gets absorbed by 
the other membrane. In the intermediate steps this third membrane crosses the branch disk 
from the first R^ into the second, flipping its orientation. The self-intersection involved in 
the flipping of the membrane orientation is shown in hgure [T^ which is an enlargement of 
the small rectangle marked in the fourth panel of figure [TH 

Figure [T^ depicts a case intermediate between the previous two examples, corresponding 
to 9q = —7r/6. In this case the transfer of angular momentum does not involve the exchange of 
a third membrane, but the interplay between the equipotential surfaces is rather interesting. 
In figure [T4| we show a detailed view of the area within the rectangles marked in the third and 
fourth panels in figure 13, presenting a sequence of contour plots for values of the potential 
between cp = 0.380 and cf) = 0.378. The fourth panel in figure 14 shows clearly that in this 


case the splitting takes place not at a point, but simultaneously at points along a ring. This is 
of course not a generic feature. It is a degenerate case which only occurs in axially symmetric 
solutions for certain values of the parameters. A similar splitting ring is also present in the 
previously discussed case of hgure 

It is quite remarkable that the rather striking behaviour of membranes illustrated by the 
above examples may be reproduced by solutions of such a simple and well studied equation 
as the Laplace equation. 


4.3.2 Comments on the moduli space of solutions 

In this subsection we discuss some features of the moduli space of instantons associated with 
general Riemman spaces. As we previously noted, the number of copies of R^, re, corresponds 
to the number of membranes in the hnal state. Positive charges, Ji > 0 {i = 1,..., re), can 
only be placed at the origin of these spaces. The charge Ji corresponds to the angular 
momentum of the z-th membrane in the initial state. Some of the Jj’s may be zero and the 
number of non-zero JJs is the number of membranes in the initial state. 
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Figure 9: Evolution of the membrane profiles for a process with two membranes in the initial 
state and two in the final state {Oq = —27r/3). 
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Figure 10: Enlargement of the area marked by a square in 


the third panel of figure 


The sheets of the Riemann space are connected by branch disks, bounded by branch loops. 
By definition when a branch disk connects two copies of IR^ the disk is located at the same 
place, with the same shape, in the two copies Except for this condition, the number of 


branch loops connecting the copies of IR^, their positions and their shapes are arbitrary. In 
the following we collectively refer to the number of copies of IR^ and the number, the positions 
and the shapes of the branch loops as the geometry of the Riemann space. Our conjecture is 
that for any geometry of the Riemann space (and the distribution of positive charges at the 
origins) we have a unique solution to the Laplace equation satisfying the boundary conditions. 

A point on the instanton moduli space is characterised by specifying the geometry of the 
Riemann space. This will give an interpretation of the moduli space of solutions of the BPS 
instanton equation (2.12) discussed in 18 when J is large (but finite). Stated differently. 


the moduli space associated with (2.12) can be considered as a regularisation of the moduli 
space of branch loops. 

The outgoing flux at infinity in the i-th copy of IR^, which we call J', corresponds to the 
angular momentum of the f-th membrane in the hnal state. 

A necessary and sufficient condition for the existence of solutions to the instanton equa¬ 
tions interpolating between the initial state characterised by Jj, i = 1,... ,n, and the final 
state characterised by J', i = l,...,n, was established in 


18 


Below we will show that 
the same condition can be derived from our approach taking advantage, in particular, of the 
linearity of the Laplace equation. 

The linearity implies that for a fixed geometry of the Riemann space, one has a linear 
relationship between Jj and J', 

j[ = Y,K^Jj. (4.23) 

j 

The matrix K characterises the flow of the flux lines for a given Riemann space. (It is 
reminiscent of the coefficients of capacity in the electrostatic theory of conductors, which are 
determined by the shape and positions of the conductors.) In general the matrix K is not 
symmetric. We will also use a matrix-vector notation 


J' = KJ. 


(4.24) 


'^If this condition were not met, we would get unphysical discontinuities in the instanton solutions in the 


j/-space, reconstructed using (3.7l. 
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Figure 11: Evolution of the membrane profiles for a process with two membranes in the initial 
state and two in the final state {Oq = —vr/lO). 
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Figure 12: Enlargement of the area marked by a rectangle in the fourth panel of figure |ll| 
where J = (Ji, ..., J^), J' = (J{,..., J^). 

Some simple but important properties follow from the definition of Kij. The first property 
is the positivity of the elements of the matrix K 


Kij > 0 . 


(4.25) 


In order to see this, we consider the case in which there is only one unit charge, located at 
the origin of the j-th copy of IR^, in the entire Riemann space. The flux flowing to the i-th 
space from the j-th space equals K^. Since there is a positive charge in the j-th space and 
no charge in the i-th space, it is clear that the flux should always flow from the j-th to the 
i-th space, not the other way around. This implies Kij > 0. 

Other important properties of K are the sum rules. 


= l Vj, 


(4.26) 




= 1 y i. 


(4.27) 


The property (4.26) simply follows from the conservation of the number of flux lines (i.e. 
the Laplace equation and Gauss’ theorem). The property (4.27) can be deduced as follows. 


We consider a special potential function in the Riemann space defined by the requirement 
that in each copy of IR^ it is identical to the Coulomb potential associated with a unit charge 
at the origin. Such a potential solves the Laplace equation and satishes all the boundary 
conditions The existence of this solution implies that if Ji = (1,1,1,..., 1), then J[ = 
(1,1,1,..., 1) in (4.23), which is equivalent to (4.27). 


The positivity (4.25) and the sum rules (4.26)-(4.27) imply rather strong constraints on 


the allowed combinations of Jj’s and J-’s, including conditions equivalent to the criteria given 


m 


18 for the existence of solutions. 


First, one can show that 




J‘ 


IP 


<T,J! 


(4.28) 


We note that this is a generalisation of the observation given below (4.211. 
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Figure 13: Evolution of the membrane profiles for a process with two membranes in the initial 
state and two in the final state {Oq = —vr/O). 
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Figure 14: Enlargement of the area marked by rectangles in the third and fourth panels of 
figure 13 showing intermediate steps between the two. 
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for any p > 1. By using the vector norm ||n||p = the matrix norm is defined 


by 32 


= max- 


v^o ||n| 


Hence in order to establish (4.28) it is sufficient to show 

ll^llp< 1- 

Using an inequality for matrix norms (cf. formula (1-11) in |33|), 

and the well-known formulae 1321 


1 = max > IA 
i=l 




= max 

l<2<n 




(4.29) 


(4.30) 


(4.31) 


(4.32) 


j=i 


we see that Ijii'lli = 1 and ||ill||oo = 1 and hence ( |4.30[ ) and (4.28) are established. For p = 3, 
we have 

(4.33) 


which gives a continuum counterpart of the necessary condition (2.16) for the existence of 
instantons (which can also be derived directly by defining the continuum version of (2.11)). 

One can introduce an “entropy” associated with the distribution of the total angular 
momentum among the individual membranes defined as 


^ (-Ji log Jj) 


and 




(4.34) 


for the initial and final states respectively. Then by taking p = 1 -|- e in (4.28) and using the 
conservation law = X] *4^ one can show that the entropy is always non-decreasing, 

J](-J'logJ') > J](-jaogJ.). (4.35) 

i i 

Let us now deduce the necessary and sufficient condition for the existence of instantons given 
in 18 . The condition can be phrased as follows (cf. figure 3 in [18]). Since the Jj’s are 


non negative, they can be represented by a histogram as shown in figure 15 (i). A similar 
histogram can be constructed with the JLs in the final state. We draw a horizontal line at a 
height J > 0 on the histogram and denote the area of the histogram below the value J by 


A(J; J). 


The condition in 18 can then be expressed as 

A{J]J')>A{J]J) VJ. 


(4.36) 


(4.37) 


(We note that we do not assume Ji > J 2 Jn or > J 2 >...> Jn here.) Thus, 

the area of the support of the histogram is always larger than or equal to that of the 
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Figure 15: The angular momenta of membranes J = (Ji,..., Jn) can be represented as a 
histogram as shown in (i). The area of the histogram J below a given height J is denoted by 
A{J]J). We define histograms J(a) and J(b) by cutting the original histogram J along the 
horizontal line and taking the part above and below it, as shown in {ii) and {in). 


histogram Ji. An interpretation of this condition is that the histogram J' should be more 


smeared compared to the histogram Jj, which is natural in view of (4.35) 


12 


In order to prove (4.37) we introduce two auxiliary histograms, and generated 

we define J(a)i and 


15 


cutting the histogram Ji at a fixed height J. As shown in figure 

as the histograms above and below the horizontal line at the height J, respectively. By 
definition we have 


J — "^{a) T 


l{b)- 


(4.38) 


Moreover 

J' = KJ = KJ^a) + KJ(^b)- (4.39) 

We first focus on the histogram defined by KJ(by All columns in this histogram have height 
smaller than or equal to J, i.e. {KJ(^b))i ^ J- To see this we recall that (4.27) implies that 


K 


■ J' 


■ J ■ 

_ J 


J _ 


(4.40) 


Since all elements of Jy,-^ are smaller than or equal to J by definition and all elements of K 
are positive (4.25), the elements of KJ^p-^ cannot be larger than J. 

The sum rule (4.26) implies that area of the histogram KJ^b) is the same as the area of 


the histogram J^b)^ which by definition equals M(J; J). Now recall that J' = KJ(^^-^ + KJ(py 
Since the height of the histogram KJ(^b) is less than or equal to J, all of the histogram KJ(^b) 
contributes to A{J', J'). The area A{J', J') has a further contribution from KJ(^a)^ which is 
always non-negative, because of the positivity of K and J[a)- Thus, we have 


A{J-,J')>A{J;KJ^b))=A{J-,J), 


(4.41) 


and (4.37) is proven. 


We note that the condition (4.371 is stronger than the condition (4.331. For example, the case with 


J = (5,5,0) and J' = (6,2,2) is allowed by (4.331, but it is prohibited by (4.371. 
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Another interesting property of the moduli space of solutions to the BPS instanton equa¬ 
tion (2.12), which was established in [18] , is the additivity rule (2.24) for the dimension of 
the moduli spaces. A possible interpretation of this result, which is instructive if perhaps not 
mathematically rigorous, is the following. We assume that we have two instanton solutions, 
interpolating between vacua A and B, and between vacua B and C, such that both of them 
are reasonably localised in time. It is natural to expect that one can construct an instanton 
solution interpolating between vacuum A and vacuum C, by stitching together the two origi¬ 
nal solutions taken to be well separated in time. The number of parameters of the composite 
instanton thus constructed should be given by the sum of the numbers of parameters of the 
two original instantons. 

This construction of a composite instanton from two instantons well separated in time has 
a natural counterpart in our description of the instantons using Riemann spaces. The follow¬ 
ing argument relies on viewing a solution in a multi-sheeted Riemann space as a composite 
of solutions built using pairs of connected copies of IR^. Within this picture the basic idea 
underlying the construction is that, at least in the asymptotic region at infinity, a branch disk 
with incoming flux can be viewed as a point charge. Here we wish to point out the essential 
features of this construction without going into the details. In the following it is crucial that 
the time variable in the BPS instanton corresponds to the logarithm of the distance from the 
origin in our description. We start with two instanton solutions, each associated with the 
branch disks Di and 02 - Because of the scale invariance of the Laplace equation, one can 
move one of the disks, e.g. D 2 , further away from the origin by rescaling simultaneously its 
size and its distance from the origin. The physics associated with the second disk, such as 
the number of flux lines going into the disk and the profile of the potential function do not 
change under the rescaling. We refer to the rescaled disk as Now we consider a new 
Riemann space equipped with two branch disks, Di and Since the disk D 2 is sufficiently 
far away, its existence does not affect the potential produced by the point charge at the origin 
and the first branch disk Di. Furthermore, from the viewpoint of the second disk D 2 , the 
disk Di and the point charge appear as a new point charge at the origin whose value can be 
reconstructed using the original charge and the number of flux lines going through the disk 
Di. In this manner, the potential on the Riemann space with two disks {Di, D' 2 ) can be 
constructed by using the potentials corresponding to the two separate solutions associated 
with the disks Di and D 2 . 

As previously mentioned, the calculation of physical transition amplitudes using a semi- 
classical approximation involves the integration over the collective coordinates parameterising 
the instanton moduli space. From the point of view of the description of the moduli space in 
terms of the geometry of Riemann spaces, this means that one should sum the contributions 
of all possible combinations of branch surfaces giving rise to the same initial and final states. 
Such a sum may include not only branch disks of more complicated shapes than the circular 
ones we discussed, but also in principle branch surfaces with more complicated topology, e.g. 
consisting of higher genus surfaces with boundaries. It would be interesting to study the 
characteristics of the resulting Riemann spaces. 
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5 Conclusions and discussion 


In this paper we have initiated a study of the interactions associated with membrane split¬ 
ting/joining processes in M-theory. We considered a large angular momentum sector of M- 
theory in AdS 4 x 5^, which can be described using the matrix model approach in the pp-wave 
approximation. The vacua of the pp-wave matrix model consist of collections of concentric 
membranes and contributions to splitting or joining processes arise from tunnelling ampli¬ 
tudes between states built on vacua with different numbers of membranes. Classical tunnelling 
configurations are instantons, i.e. finite action solutions to the Euclidean equations of mo¬ 
tion of the pp-wave matrix model. Relying on the same ideas underlying the definition of the 
matrix model as regularisation of the supermembrane theory, we have presented a method 
for the construction of approximate solutions to the instanton equations. Using a continuum 
formulation in which large matrices are approximated by functions, we have shown that the 
problem of constructing solutions to the instanton equations of the pp-wave matrix model 
can be reduced to that of solving the three-dimensional Laplace equation. A remarkable 
outcome of our analysis is that in order to describe membrane splitting/joining processes it 
is necessary to study the Laplace equation in a Riemann space, which is a three-dimensional 
generalisation of the familiar concept of Riemann surface. The interpolating configurations 
at different points in Euclidean time arise as equipotential surfaces for the solutions to the 
Laplace equations in a Riemann space with suitable boundary conditions. The use of Rie¬ 
mann spaces provides a very compelling answer to various issues arising in the description of 
membrane splitting processes in terms of the Laplace equation. 

The approach we developed in this paper presents some similarities to the description of 
splitting/joining interactions in the matrix string formulation of type IIA string theory 
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in terms of instantons in two-dimensional SYM theory 35f(37 
between our work and k 


There are also some analogies 


, in which 1/4 BPS solitons in AA = 4 SYM are described using 
multi-pronged {p, q) string states. The latter are constructed using the two-dimensional 
Laplace equation which arises from configurations of M2-branes stretched between M5-branes. 

In section we presented exact solutions to the Laplace equation, which represent ap¬ 
proximate solutions to the original BPS instanton equations valid for large J. The instanton 
configurations in the matrix model are saddle points for the Euclidean path integral of the 
theory and can be used to compute contributions to physical transition amplitudes using a 
semi-classical approximation. As a next step it will be important to deduce the exact form 
of the solutions in the matrix model, associated with the approximate solutions we obtained 
using the continuum approximation. The connection between the instanton equations and 
the Laplace equation in the continuum approximation is established mapping the former to 
the Nahm equations. This means that in the process of reconstructing the matrix model 
solutions we should be able to take advantage of the integrability properties of the Nahm 
equations. 

Semi-classical contributions to physical transition amplitudes in the Lorentzian signature 
pp-wave matrix model should be obtained from a saddle point approximation about the 
instanton configurations followed by an analytic continuation. The semi-classical calculation 
involves an integration over the moduli space of the BPS instantons, whose properties were 
studied in EZlI Ts] . Completing the semi-classical analysis will allow us to extract effective 
interaction vertices for membranes in the pp-wave background. A first interesting issue 
to address is the determination of the associated coupling constants. Naively one would 
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expect instanton induced transition amplitudes to be exponentially suppressed, due to the 
contribution of the classical instanton action. For instance the most elementary amplitude, 
corresponding to the splitting of a single membrane with angular momentum J into two 
membranes with angular momenta Ji and J 2 , carries a weight 


c _ JJ-1J2 

^ OE — Q 


(5.1) 


where Ji ~ J 2 ~ J and <€. J <€. for the applicability of the pp-wave approxi¬ 

mation in a weak coupling regime. However, it is probably premature to conclude that the 
coupling constant measuring the strength of membrane interactions is exponentially small. 
The dependence on the parameters can be modified by factors arising from the integration 
over the instanton moduli space, which has a dimension that grows with J - for example it 
is 4 J 2 in the case of the single membrane splitting, see (2.23). Volume factors can modify 


the exponential behaviour (5.1) and it is possible, at large J in particular, for the complete 


physical transition amplitudes to have a power-law dependence on the parameters J and N, 
which might be more natural from the point of view of a description of membrane inter¬ 
actions in target space. A similar mechanism, leading to a power-law behaviour, has been 
recently suggested to arise in the calculation of solitonic contributions to certain scattering 
amplitudes 39 and may be relevant for the understanding of the connections between the 


maximally supersymmetric Yang-Mills theory in five dimensions and the (2, 0) superconfor- 


mal theory in six dimensions 40 41 


The particular solution (2.21) constructed in 18 


which approaches the configuration 
corresponding to the J-dimensional representation l©l©---©latt = + 00 , was studied in 
the context of the pp-wave matrix model in 17 . According to a proposal presented in |42|, 


the configuration associated with the representation 1 © 1 © • • • © 1 should be interpreted 
as a single M5-brane and therefore the solution (2.21) would contribute to an amplitude 


coupling M2- and M5-branes. It would be interesting to determine whether such a coupling 
has a different weight in terms of the J and N parameters, compared to those involving only 
M2-branes. 

The study of instanton effects in another maximally supersymmetric theory, AA = 4 SYM 
with SU(V) gauge group, suggests that the large J limit, which plays a central role in our 
continuum approximation, may also lead to important simplihcations in the semi-classical 
calculation of physical transition amplitudes. In the case of AA = 4 SYM at large N the 
integration over the multi-instanton moduli space is dominated by special configurations and 
a saddle point approximation makes it possible to compute the leading contributions at large 
N for arbitrary instanton number |43| . It would be very interesting to understand if a similar 
saddle point based approach can be applied in the present case for large J. 

The following considerations, albeit somewhat speculative, may provide insights into pos¬ 
sible mechanisms leading to the emergence of a saddle point approximation for large J. 
Among the collective coordinates parameterising the instanton moduli space there are vari¬ 
ables associated with the SU(J) gauge orientations, which are the remnant in the matrix 
model of the invariance of the membrane theory under area preserving diffeomorphisms 
(APD’s). When computing contributions to gauge invariant observables the integrations 
over these colour-related variables can be carried out and in general produce a non-trivial 
measure for the integrals over the remaining collective coordinates, that can be referred to 
as the ‘physical’ ones. At large J this measure can be very peaked and give a large weight 


41 















to special regions of the physical moduli space. In such a situation it is possible to expect 
the emergence of a good saddle point approximation. Given the fact that the SU(J) gauge 
symmetry descends from the APD invariance in the continuum, it may be natural to expect 
that the configurations acquiring the largest weight in the physical moduli space should be 
the smoothest ones, since they are the ones for which the hnite SU(J) symmetry is closer to 
the infinite dimensional group of APD transformations. In our description of instanton con¬ 
figurations in terms of the Laplace equation on a Riemann space the collective coordinates on 
the instanton moduli space are associated with the geometry of the Riemann space. Then the 
above reasoning suggests that contributions arising from the smoothest branch loops may be 
the dominant ones. Hence for the specific transition amplitude corresponding to the splitting 
of one membrane into two, the Riemann space with circular branch disk used in the Hobson 
solution that we employed in section [^should presumably provide the dominant contribution. 

Supersymmetry of the pp-wave matrix model is important in various aspects of the cal¬ 
culation of instanton contributions to the transition amplitudes that we are discussing. It is 
expected to be responsible for the cancellation of the determinants arising from the integra¬ 
tion over non-zero mode fluctuations in the instanton background. Moreover the instanton 
moduli space includes fermionic collective coordinates associated with supersymmetries bro¬ 
ken by the instanton configurations. Although the particular features of the supersymmetry 
algebra in the pp-wave background make the analysis rather subtle, we expect the presence 
of fermion zero modes to restrict the set of allowed instanton amplitudes. In addition to the 
conditions for the existence of interpolating instanton configurations discussed in |18| , further 
selection rules should arise from the integration over the fermionic collective coordinates. This 


was shown in |17| to be the case for the special solution (2.21) obtained in |18| . In general the 
results of 


showing that the vacua of the pp-wave matrix model are non-perturbatively 
protected, indicate that instanton induced transition amplitudes can be non-zero only when 
excited states are involved. It should be possible to gain further insights into the general 
selection rules which constrain these processes by studying the fermion zero modes in the 
background of the general matrix model solutions corresponding to those discussed in this 
paper in the continuum approximation. Similarly to what happens in the case of AA = 4 
SYM it should be possible to carry out this analysis even without a complete calculation 
of the transition amplitudes. 

Various properties of the moduli space of solutions of the BPS instanton equations were 
discussed in dill T^ . The collective coordinates parameterising the matrix model solutions 
have counterparts in the parameters charactering the solutions to the Laplace equation that 
we used our approach. The description in terms of the Laplace equation provides a ge¬ 
ometric and intuitive interpretation of some of these parameters and the linearity of the 
equation is useful in explaining certain features of the solutions. According to our proposal 
the configurations interpolating between different vacua of the pp-wave matrix model have 
an approximate description in terms of equipotential surfaces of potential functions solving 
the Laplace equation in certain Riemann spaces. Any given solution is characterised by the 
states at t = Too in the original picture, i.e. by the asymptotic SU(2) representations de¬ 
fined by the dimensions, Ji, i = 1,..., m, and J^, /c = 1,..., n, of the irreducible blocks at 
t = —oo and t = -(-oo respectively. In our description the values of the Ji integers for the 
initial state correspond to point charges at the origin in different copies of while the J^’s 
for the hnal state correspond to the fluxes of the gradient of the potential at infinity in each 
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of the R^’s which constitute the Riemann space. In such a picture the information on the 
parameters characterising the solutions should be encoded in the geometry of the Riemann 
space, e.g. the shape and relative position of the branch disks. In section ^ we discussed 
some properties of the moduli space of the BPS instantons from the point of view of the 
Riemann space picture. It would be interesting to have a more complete understanding of 
the results of |18| within the framework proposed in this paper. In |18| , it was shown that 
the BPS equations, with the same boundary conditions considered in this paper, describe 
domain wall solutions in the N = 1* SYM theory and also the D3-D1 system in the context 
of the AdS/CFT correspondence. It would be interesting to consider the implications of our 
results in these contexts. 

The matrix model, providing a regularisation of the supermembrane theory, represents 
the proper framework for a full quantum mechanical calculation of physical transition am¬ 
plitudes of the type we discussed in this paper. However, it may be possible to carry out 
the semi-classical calculation using as starting point the path integral for the continuum su¬ 
permembrane theory. In the semi-classical approximation the calculation is reduced to an 
integration over the instanton moduli space. It may then be possible to introduce an appro¬ 
priate regularisation, suitable to make the integration over the moduli space of the continuum 
solutions well defined. In this way an approximate result, valid for large J, might be attained 
working directly with our continuum solutions without reconstructing the explicit form of the 
corresponding matrix model configurations. Such a possibility is made more plausible by the 
large amount of supersymmetry in the theory. The use of this continuum approximation can 
probably be justihed more straightforwardly if indeed the integration over the moduli space 
is dominated by a saddle point for large J. 

The study of interactions between M-theory objects using the matrix model formulation 
has so far been mostly limited to processes involving no longitudinal momentum transfer. 
These calculations allow to extract effective supergravity interactions from perturbative cal¬ 
culations in the matrix model. A comprehensive review of these results and a detailed list of 
references can be found in [45] . A contribution to the scattering of membranes with minimal 
momentum transfer in the M-theory direction was carried out in 46 . The process studied 
in 46 - described in terms of transfer of a DO-brane between two parallel D2-branes - can 


be considered as analogous (in a special kinematical regime) to the one corresponding to 
the solution with two membranes in the initial state and two in the hnal state discussed in 
section [4^ - more specifically the example shown in figure |TH However, we emphasise that 
the approach proposed in this paper is more general, since it allows us to describe processes 
with an arbitrary number of membranes in the initial and final states. Moreover we can 
describe interactions involving the exchange of a large amount of longitudinal momentum. In 
particular, in cases such as the example shown in figure]^ we have an exchange of a genuine 
M-theory object, which in the DO-brane picture would require considering a configuration 
involving a large number of DO-branes. In a calculation of the type presented in |46j this 
would mean evaluating contributions from sectors with large instanton number. Further work 
on interactions involving longitudinal momentum transfer was done in [47] . 

Our results on membrane splitting/joining transitions provide a foundation for the sys¬ 
tematic study of more general interactions of membranes in M-theory. The existence of 
classical solutions corresponding to configurations that interpolate between states containing 
different numbers of membranes, supports the conclusion that the matrix model is capable 
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of describing splitting/joining interactions, reinforcing the early indications offered by the re¬ 
sults of 46,47 . To provide a truly convincing argument it is of course necessary to complete 


the calculation of the physical transition amplitudes in semi-classical approximation. More¬ 
over, it is important to be able to check the results by independent means and the AdS/CFT 
duality of provides a way of achieving this. Based on the dictionary for the large J sector 
established in |^, the vertices associated with membrane splitting/joining interactions can 
be related to correlation functions of monopole operators in the dual CFT. As pointed out 
in [^, the use of radial quantisation is a convenient tool for the study of monopole operators 
and their correlators in the ABJM theory. The relevant monopole operators are analogous to 
disorder operators in two-dimensional CFT. To compute CFT quantities that can be related 
to membrane interaction vertices one needs to consider amplitudes involving multiple inser¬ 
tions of such operators. It is an interesting possibility that Riemann spaces may be useful 
for the analysis of correlation functions of this type in the radially quantised ABJM theory. 
Moreover, as is the case for the study of the spectrum [^, the presence of the large parameter 
J should make it possible to develop an approximation scheme allowing a reliable comparison 
with the matrix model results. A similar approximation, valid in a sector characterised by a 


large global charge J, was recently discussed in 48 in the case of a different class of strongly 
coupled three-dimensional theories. 

An alternative way of comparing the effective membrane interaction vertices to the dual 


CFT may be based on the method proposed in 49 for string theory in a pp-wave background. 


In this approach membrane couplings would be related to matrix elements of the dilation 
operator of the ABJM theory between states associated with monopole operators. It should 
be possible to compute such matrix elements in semi-classical approximation, using tunnelling 
configurations interpolating between monopole states characterised by different GNO charges. 
It is natural to expect that the Nahm equations might play a role in a gauge theory calculation 
of this type as well. We note that certain BPS equations associated with the description of 
bound states of M2- and M5-branes in the ABJM theory were studied in 50,51 , where 
interesting relations to the Nahm equations were pointed out. 

It would be very interesting to generalise the type of analysis that we developed for the 
pp-wave matrix model to other cases, relevant for M-theory in different backgrounds. The 
cases of the AdS 4 x 5’^ and AdSy x spaces are of course of interest. For these backgrounds 
no matrix model formulation is known, however, the supermembrane action has been con¬ 
structed [^. Therefore at least a study of instanton solutions using our approach might 
be possible, since the continuum version of the instanton equations can be also understood 
directly from the membrane theory without reference to the matrix model as discussed in 
section 3.2 Another obvious case to consider is flat space, which may be approached 
starting from the general pp-wave matrix model described in appendix and using the mass 
parameter ^ as an infra-red regulator. The /r —)• 0 limit would, however, be very subtle and 
require extreme caution. 

The emergence of the notion of Riemann space in connection with membrane interactions 
is very suggestive. We were led to introduce this concept in order to describe the evolution 
of membrane configurations in terms of solutions to the Laplace equation. It is an intriguing 


For the membrane instantons in the AdS space, solutions with run-away behaviour (discussed at the end 


0 


of section 4.1 ) may be relevant in view of the interpretation of the change of variables mentioned in footnote 


44 















possibility that our results may be an indication of more general features of the dynamics 
of membranes and that Riemann spaces may turn out to be a central ingredient in the 
description of membrane interactions, in the same way as Riemann surfaces are essential in 
the formulation of string perturbation theory as a genus expansion. 
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A Geometric proof of the equivalence between Laplace and 
continuum Nahm equations 


In this appendix we provide a proof of the equivalence between the three-dimensional Laplace 


equation and the continuum Nahm equations (3.8). The latter are sometimes also referred 


to as SU((X)) Nahm equations. The argument that we present is more geometric than those 
in |20l|21|. We start with the continuum Nahm equations, 




(A.1) 


where a is a positive constant and the Lie bracket is defined as {f,g} = ~ ^r^^ 

for arbitrary functions, / and g. This equation describes the evolution in time s of a two- 
dimensional surface parametrised by (ci^, cr^). Locally this evolution defines a function 4>{z^) = 
s{z^) on a region of the three-dimensional 2 ;-space swept by the evolving surface. We prove 
that the function (f> = s satisfies the Laplace equation. 

It is useful to rewrite the continuum Nahm equations in the form 


dz 

ds 


d 


d 




(A.2) 


This can further be rewritten as 

dz 
~ '¥s 


d 


d 


dS 


“dcTldu2- 


(A.3) 
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Here we considered an infinitesimal area element dcr^dcr^ around cr^) in cr-space and the 
corresponding area element in 2 ;-space, dS, and used the fact that 

From the definition of the gradient it follows that 

ds = d(/> = V(/> • dz, (A.5) 


and hence 




(A.6) 


Comparing (A.3) and (A.6) we obtain 


V(j)-dS= -duMa^. 
a 


(A.7) 


Applying this formula to the volume element depicted in figure 16 which is a tube whose 
sides are constructed from the electric flux lines associated with cj) and the bottom and top 
plates are the area elements corresponding to dcr^dcr^ at (/) = s and (/> = s + ds, we see that 


j V(p-dS = 0 (A.8) 

since the sides of the tube, where • dS^ = 0, do not contribute and the contributions of 
the top and bottom cancel each other. This implies the three-dimensional Laplace equation. 


One can also derive the continuum Nahm equation from the Laplace equation. Starting 
from 4>{z) satisfying the Laplace equation, we construct z{s, cr^, by solving (j){z) = s. This 
requires fixing the cr^, ci^ coordinates. We first fix cr^, for a given s by the requirement that 
the flux density be constant as in (A.7). The definition of cr^, ci^ can be extended for arbitrary 
s, by requiring that ^ be proportional to V(/). Once the paramterisation is fixed, one can 
deduce the continuum Nahm equations from the Laplace equation by basically backtracking 
the above proof. 


B Behaviour at the splitting point 


In this appendix, we discuss the behaviour of families of equipotential surfaces for solutions 
of the Laplace equation of the type presented in section near points where the splitting 
occurs. 

The splitting point is a singularity of the continuum Nahm equation, which is a partial 
differential equation (PDF). It is a well-known idea to study such singularities by describing 


the solutions of a PDF as a level set of a given function 53 . The singular point of the PDF is 


not singular for the function and thus this approach provides a way to discuss the phenomena 
described by the PDF avoiding the singularity. In the present case the PDF describes the 
evolution of a surface and another representation of the same membrane evolution, also free 
of singularities, is provided by the matrix regularisation. The equivalence between these two 
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Figure 16: An infinitesimal flux tube. The area element els' defined in (A.4) is proportional 


dz 


to both V(/> and 


regularisations - the level set description of the PDF and the matrix model formulation - is 
implicitly assumed in our analysis. 

In a region sufficiently close to a generic splitting point, the potential cj) can be approxi¬ 
mated by 

(B.l) 


s = cj){x, y, z) = -{x"^ + y2) _ 2 :^, 
which satisfies the Laplace equation, up to suitable rescaling of the x, y, z variables. The 


surface described by (B.l) is a connected hyperboloid, a cone or a disconnected hyperboloid. 


depending on whether s>0,s = 0ors<0 respectively. The s = 0 surface embedded in the 
(x, y, z) space is a cone with tip at z = 0. The continuum Nahm equations are 


X = -a{y, z}, ij = -a{z, x}, i = -a{x, y}, 


(B.2) 


where a > 0 and the Lie bracket is defined in (3.4). At the tip of the cone the equations 


become singular, while the potential (B.l) itself is, of course, non-singular 


The form of the potential (B.l) can be used to study universal features of the splitting 


point. In particular, in order to discuss the relation between the matrix regularisation and 
the regularisation via the Laplace equation, one can assume the matrix size to be sufficiently 
large. Hence one can focus on the vicinity of the splitting point, i.e. the tip of the cone, and 


the generic form of the potential (B.l) is sufficient. 


We discuss the relation of the two regularisations in the following way. Firstly, we rewrite 


the solution to the Laplace equation near the singularity (B.l) using a hodograph transforma¬ 


tion, i.e. exchanging the role of independent and dependent variables, to obtain x(s, cr^, cr^). 


*In this appendix we parametrise what was referred to as the z-space in the main text by {x, y, z). 
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y{s, cr^, fj^) and z(s, a^,(T^). As functions of the (s, cr^, cr^) these satisfy the continuum Nahm 
equation (B.2), see [^[M] and appendix [A| It turns out that this solution is characterised by 
a cubic equation. Then we discretise the functions, x(cr^,cr^), y(cr^,iT^) and z((T^,(T^), using 
the matrix regularisation and check if these satisfy the Nahm equations when the matrix size 
is large enough. 


To construct x, y and 2 ; satisfying (B.l) as functions of we proceed as follows. 


For convenience we introduce a radial coordinate p, 

p = + 7/2 _ x/2x/^2 + s . 


(B.3) 


where in the last equality we have used (B.l). We fix the a coordinates by using ( A.7| ). We 
compute the infinitesimal area in cj-space as follows. Using 

d(Area in (x, y, 2 ;)-space) = 27rp^ 1 + — 2vr-\/2-\/32;2 + s dz (B.4) 


and 


I V(/(>| = |(x, 7/, —22;)! = -\/x2 + 7/2 _|_ 42;2 = 0^2 + 2s = '/2\/3z'^ + s, (B.5) 

we obtain 

d(Area in cr-space) = Q:47r(32;^ + s)d 2 ;. (B.6) 

The infinitesimal area above refers to the portion of the membrane between 2 ; and z dz. 

We parametrise the cr-space by ^ and ip with range (— 00 , 00 ) and [0, 27r) respectively. We 
have chosen the parametrisation such that ^ increases as 2 ; does and (p is the standard polar 
angle. 


Integrating both sides of (B.6) we get a cubic equation, 

q;(2;^ 


e 




(B.7) 


from which we obtain 2 ; (and p via (B.3)) as a function of ^ and s. When there exist three 


solutions in (B.7), we choose a branch: for positive (negative) a positive (negative) root 
of the equation is employed. For s > 0 there is only one root for every ^ G (— 00 , 00 ). For 
s < 0 there is a discontinuity at ^ = 0, such that 2 ; = and 2 ; = for ^ —)• 0 from 

above and ^ —)• 0 from below, respectively. Therefore, the cr-space is cylindrical for s > 0 
(corresponding to a connected hyperboloid), whereas it consists of two infinite planes for 
s < 0 (corresponding to a disconnected hyperboloid) and a cone for s = 0. It is interesting 
that a family of hyperboloids (described for each s by a quadratic equation) is parametrised 
in this way by a cubic equation. 

The ip coordinates thus defined can be identified as the (cr^, cr2) coordinates directly so 


that we have 


and 


{/>5} = 


Of dg 


dip 


dcjMcj^ = 


dip 

d^d(/7. 


(B.8) 


(B.9) 


One can study the singularity at the splitting point, using the coordinates ^,ip. A natural 
measure for the singularity is the “energy density” 


^ = \ ({a;, y}^ + {y, + {z, x}^) 


(B.IO) 
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We use the relation between the continuum Nahm equation and the Laplace equation to 
obtain 

1 1 1 

(B.ll) 


e = 


a2|V(/>P + 4^2) (6z‘^ + s) 

Thus the “energy density” diverges when s = 0 at z = 0 (i.e. ^ = 0) as 


(B.12) 


We note however that the divergence is weak enough so that the integral of the “energy 
density” 

j ed^d^ (B.13) 

is finite. 


The coordinates, x, y and z, are functions of (s,^, (/?); from (B.7) z{s,^,(p) = z{s,^). The 
variables x and y can be conveniently described by using 


w = X + iy = = x — iy = pe 


(B.14) 


It is easy to check that the functions z and w satisfy the continuum Nahm equations (B.2), 


z = —-a{w,w^}, w = ia{w,z}, 


(B.15) 


where we have introduced i = dgZ and z' = d^z. The above equations are equivalent to 

(B.16) 


i = —ap'p, p = apz'. 


Acting with dg and on both sides of the cubic equation (B.7) we have 

, 1 


z = — 


z' = 


+ s ’ 2 a{2>z‘^ + s) 


Using (B.3), these imply 


P = 




y/2(3z^ + s)’ ^ y/2a\/ z"^ + s{3z^ + s) 


(B.17) 


(B.18) 


Substituting into the continuum Nahm equations (B.16) we see that the equations are satis¬ 
fied. 

We then discretise the functions, z and w, via the matrix regularisation and we check if 
they satisfy the Nahm equations when the matrix size is large enough. Here we write the 
Nahm equations in the form 


X = -a^[Y,Z] Y = -a^[Z,X], Z =-a^[X,Y], 


or equivalently, 

where 


Z = -1^[W,W\ W = -[W,Z], 


W = X + iY, =X-iY. 


(B.19) 

(B.20) 

(B.21) 
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Here X, Y, Z and W are infinite dimensional matrices (corresponding to the non-compact 
cj-space). If the original matrix size J is sufficiently large, we have C <C 1. 


15 


Our ansatz for the matrices is as follows 27,54]. We take Z to be diagonal and assume 


that W has only non-zero sub-diagonal entries, i.e. the only non-zero components in the 
matrices are Zmm = Wm,m+i respectively. We use a slightly unusual convention in which 
the rows and columns of infinite-size matrices are labelled by half-integer valued indices. 
The index m can take both positive values, m = 1/2,3/2, 5/2,..., and negative values, 
m = —1/2,—3/2,—5/2,... Using this ansatz, we propose a candidate for an approximate 


solution to the Nahm equations (B.20). For this purpose, we introduce an auxiliary function. 


z[m], defined as a root of the cubic equation 


C 


+ sz = — m. 
2a 


(B.22) 

In case this equation has three solutions, we choose a branch in the same way as described 


below (B.7). Note that in (B.22) we have discretised ^ introduced in (B.7) as .^ = Cm, i.e. 


the mesh of ^ is taken to be equally spaced, following the Bohr-Sommerfeld quantisation |27| 


Then the approximate solution to the Nahm equations (B.20) is 


Zmm = z[m], (B.23) 

Wm ,m+l — V2y/z[m -(- l/2]2 -I- s . (B.24) 

Note that this equation is well-defined even for W_i/ 2 ,+i /2 for s < 0; z[m] has a discontinuity 
at m = 0, but z[m\‘^ does not. Furthermore, we note that W_i/ 2 ,+i /2 = 0 here. For s < 0 
there are two disconnected portions, i.e. the two sheets of the hyperboloid. The values, 
m = —1/2 and m = 1/2, correspond to these two parts and we do not expect non-zero 
matrix elements between them. 


We show that these matrices satisfy the Nahm equations (B.20) up to terms which are 
sub-leading for C <C 1. Indeed, 

[W,W^]mm = Wi^^^,-Wl_,^^ 


2{z[m + 1/2] — z[m — 1/2]) {z [m + 1/2] + z[m — 1/2]) 

^,lrn]^z[m]^- — Zmm. 
dm a 


(B.25) 


Similarly, since 


[W, Z]m,m+1 = iz[m -k 1] - z[m]) Wm,m+1 

d . 


and 


IF, 


m,m+l 


dm a ^/2{2lz[mY + s) 


2z [m -k 1/2] i [m — 1/2] -k 1 ^ z[mY + s 


IF, 




V2{3z[m\‘^ -k s) ’ 


(B.26) 


(B.27) 


^ The constant C is analogous to h in quantum mechauics. For compact membranes, C is relat ed t o the 
size of matrices J by C = [cr]/(27rj) where [cr] = J (Yu. This can be understood for example from (3.21. By 
focussing on the splitting point, we have turned the original compact cr-space into a non-compact n-space, 
and thus the matrix size J should also be taken to infinity in such a way that C is fixed. 
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we have 


IW,Z] 


m,m+l 


c ■ 

hFm m+l • 
a 


(B.28) 


The existence of the approximate solution suggests that there exist an exact solution with 
the same qualitative behaviour. It would be also interesting to find the exact solution of the 
Nahm equations within the tridiagonal ansatz made above. In that case, it is known that the 
Nahm equations are closely related to those describing the so-called Toda lattice. 

C Other solutions 


C.l BHP solution 


The BHP solution (2.21) describes an instanton interpolating between a representation 
at t = —oo and the J-dimensional representation l©l©---©latt = +oo. We consider the 
special case in which is the irreducible representation J. This corresponds to a process 

in which a membrane with angular momentum J shrinks to a point-like object. Extrapolating 
the interpretation of the vacuum associated with a direct sum of irreducible representations 
as a collection of concentric membranes, the hnal state may be thought of as consisting of J 
membranes each carrying one unit of angular momentum. However, we stress that such an 
extrapolation goes well beyond the region of applicability of the approximations used in this 
paper. Using the proposal in 42 , the final state can also be understood as corresponding to 


a single M5-brane carrying angular momentum J. 

It is not difficult to construct a solution to the Laplace equation which approximates 
this instanton configuration. To obtain a solution with the right properties, we consider the 
Coulomb potential generated by a positive point charge J at the origin with the addition of 
a negative constant, 


ct){z) = 


J 


47r|z| 


— e 


—2to/R 


(C.l) 


where z = Using the transformation (3.7) in reverse order, one can compute the 


distance from the origin of the membrane corresponding to (C.l) in terms of the variables 


y = (2/^r,r 


The result is 


\y{i)\ = 


J 


{2'kT)R^ 1 + e2{*-io)/fl 


(C.2) 


This represents a sphere with time dependent radius, r{t) = \y{t)\, given by the right hand 
side of (|C.2 ). For large values of the potential, i.e. |z| « 0, corresponding to the far past 


{t —>■ —oo), ( C.1| )-(C.2) describe a conhguration which approaches a stable sphere of radius 
r = J/{2 ttTR‘^). On the other hand, the potential (C.l) becomes zero at a finite distance 
from the origin, |z| = J/(47re“^*°/^), in the far future {t —>■ +oo), where the radius of the 
membrane, \y\, shrinks to zero. 

The explicit form of the continuum counterpart of the BHP solution, in the case in which 
LI s is irreducible, is therefore 

(-oo) ’ 




(C.3) 
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where the membrane radius is 

~ (27rr)i?2 1 + e2(t-to)/R 

and n*((T^,(T^) is a unit vector in the radial direction. 


(C.4) 


C.2 Solution with positive and negative point charges 


In this section we present a new solution approximately describing the instanton interpolating 
between the irreducible vacuum J and the vacuum associated with the representation © 
1 © • • • © 1, where 1 appears with multiplicity (J — J'). From the membrane point of view 
in this process a single membrane with angular momentum J splits into a membrane with 
angular momentum J' and a second membrane which shrinks to a point in the hnal state. 

The corresponding solution to the Laplace equation is the sum of two Coulomb poten¬ 
tials associated respectively with a positive charge J at the origin and a negative charge of 
magnitude (J — J') at a point z = zo(/ 0), 




J 


J-J' 


(C.5) 


47r|z| 47r|z — Zq] ’ 

where J' is assumed to be positive. The evolution of the membrane profile, corresponding to 
the equipotential surfaces of (C.5|), can be described as follows. The relevant equipotential 


surfaces are those with s = (j) € (+oo, 0] which corresponds to t G (—oo, + 00 ). For t ~ — 00 , 
i.e. s = (/) ~ + 00 , we have a single spherical membrane of small radius. As t increases 
this membrane grows and deforms, until it splits into two membranes (both with spherical 
topology) at a value s > 0 of the potential. One of the membranes continues to grow as t 
increases towards + 00 , i.e. s = (()—>■ 0, and approaches the behaviour of the equipotential 
surfaces of a Coulomb solution with charge J'. The other membrane at ~ 0 is a sphere 
of hnite radius. (This is a well-known fact used for the electrostatic potential of a spherical 



describes the process in which a spherical membrane with angular momentum J at t = —00 
splits into a spherical membrane with angular momentum J' and a point-like object carrying 
angular momentum (J — J') 

A description of the above solution in terms of a suitable Riemann space leads to an inter¬ 
esting interpretation for the negative point charge. We consider a Riemann space consisting 
of m + 1 sheets, with m identical branch disks connecting the hrst sheet to the remaining m 
sheets. The insertion of the negative charge can be realised as the limit in which m is sent 
to inhnity while the size of the branch disks is sent to zero, so as to keep the total flux out 
of the first sheet fixed. 

It is clear that this solution can be generalised. One can put an arbitrary number of 
negative point charges at arbitrary locations. It is also possible to put the negative charges 
on different sheets of a Riemann space. 


It is also possible to consider the solution in which J' < 0. In this case, as t increases, the equipotential 
surfaces in the 2 -coordinates do not split (in the relevant region > 0) and instead approach a sphere of 
finite radius as t —>■ 00 . Thus this solution corresponds to an instanton interpolating between the irreducible 
representation J and the trivial representation 4 © ■ ■ • © i- 
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As mentioned earlier the !©• • •©! solution was conjectured in [42j to represent a single M5- 
brane. It is interesting to notice that combining this proposal with the above considerations, 
we are led to the conclusion that while M2-branes in our construction correspond to positive 
charges (which should be located at the origin of each sheet of a Riemann space), negative 
charges (which can be placed at arbitrary locations) are associated with M5-branes. 


D Relation between Hobson’s and Sommerfeld’s solutions 


The solution found by Sommerfeld 23 is 


(p{r,ro) = 


r - ro 


1 1 

—I-arcsin 

2 TT 


COS 




cosh ^ 


(D.l) 


where 


cosha5 = 




'so 


+ {zs - zsof 


(D.2) 


2rsr5o 

up to an overall constant. The notation used in the above formulae is as follows. We use the 
standard cylindrical coordinates {rs,fs,zs) for the point r where the z-axis coincides with 
the straight branch line. For the source point ro the cylindrical coordinates are (r 5 o, -z^so)- 

The solution is related to Hobson’s solution by an inversion tranformation, where the 
point of inversion A is taken to be a point on the circumference of the disk. It is clear 
that the circular branch loop in Hobson’s solution is mapped to the straight branch line in 
Sommerfeld’s solution. The relation between the two solutions was observed and numerically 
verified in work by Heise [30]. Below we provide an analytic proof. 


For the purpose of comparison, we present here again Hobson’s solution (4.13), (4.14) in 
an appropriately normalised form, 


Hr,ro) = 


1 


r - rol 


11 -I 

—I— arcsin 
2 TT ' 


COS 


Sr-Sho 


cosh ^ 


(D.3) 


where 


cosh Oh = cosh pn cosh pno — cos{ipH — ^ho) sinh pn sinh pno 


(D.4) 


and {ph: ^h) ai'e the peripolar coordinates defined in section [4.2.2 We note that we have 

used 


cosh = 


cosh an + 1 


(D.5) 


In this appendix we use the subscripts H and S to emphasise the distinction between the 
coordinates and the variables used in Hobson’s and Sommerfeld’s solutions. 


In order to show that the solution (D.3) maps to (D.l) under an inversion transforma¬ 


tion, we note that the prefactor l/|r — ro| in both solutions transforms under an inversion 


transformation in the well-known manner. Hence to show the equivalence between (D.3) and 


(D.l) under the inversion transformation it is sufficient to prove the following identities 


= Oh, 


as = an- 
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(D.6) 

(D.7) 





















The first identity needs only be satisfied up to a sign and constant shifts because the variables 
in the formula always appear in combinations such as cos {{ips — ipso)/‘^)- 

To prove (D.6) we first focus on the plane which contains the point of inversion A and the 
center of the branch disk and which is orthogonal to the branch disk. The dehnition of 6h is 
depicted in the first panel of figure We denote by P',B' the points corresponding to P,B 
by an inversion transformation around the point A. By definition we have APB = Oh- We 
also have P'B'A = cps, since the straight branch line passes through the plane orthogonally 
at the point B'. By elementary geometry, it is easy to show that APB = P'B'A, which is 
(D.6), using the definition of the inversion transformation. It is possible furthermore to show 
that the validity of (D.6) can be extended to the whole Indeed, we observe that the 
constant Oh surface in Hobson’s coordinate is a part of a sphere bounded by the branch loop. 
Since the surface contains the inversion point, it will be mapped to a half-plane bounded by 
the straight branch line. Thus a constant Oh surface is mapped to a constant (ps surface 
under the inversion and therefore it follows that (D.6) is valid generally. 

It hnally remains to show (D.7). As an and as do not depend on the Oh = fs variable, 
one can focus on the special case Oh = tt. The quantity an will then be dehned on a 
pair of points on the branch disk itself and we wish to show that it is equal to as for the 
corresponding pair of points on the branch half-plane. This can be established by noting that 
as and an are nothing but the geodesic distances in the classic Poincare upper half plane and 
disk models (the latter also referred to as the conformal disk model) of hyperbolic geometry, 
respectively 
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and hence (D.7) follows. 


These two models are related to each other by an inversion transformation 


It is interesting to note that the quantity an (or as), which is an important ingredient 
in the exact solution, has an interpretation in terms of some non-trivial geometry defined on 
the branch disk (or the branch half-plane). 


E Boundary conditions at branch disks and reformulation in 
terms of integral equations 

In this appendix, we formulate boundary conditions associated with branch disks and present 
a reformulation of the Laplace equation on a Riemann space in terms of a set of integral 
equations. The unknown variables in the integral equations are functions defined on the 
branch surfaces. Here we use the vector r = (x, y, z) to denote the position of a point in 
what was referred to as z-space in the main text. 

For concreteness we focus on the simplest case where there are two copies of IR^ connected 
by a single branch disk D. We put a point charge J at the point designated by a vector tq. 
The potential function will be denoted by </>(!)(?’) in the first space and by </>( 2 )(^") in the 
second space. 

As functions in IR^ both these potential functions have a discontinuity at the branch disk. 
In order to take care of this discontinuity, we introduce the following notation. We refer to 
the two sides of the branch disk as a and b. The values of the potential functions in each 
space, {i = 1, 2), on the a and b sides of the branch disk will be denoted respectively by 

4'ii)air) and ^(i)b(r). (E.l) 
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Figure 17: Riemann spaces 


In these expressions it is implicitly assumed that the vector r corresponds to a point on the 
branch disk which is approached from the a or b side. 

We introduce the normal vector orthogonal to the branch disk dehned at the point on 
the branch disk, n, with orientation such that n points away from the disk on the a side and 


towards the disk on the b side as shown in hgure 17 


The precise formulation of the condition that the first and the second spaces are connected 
by the branch disk is given by the following formulae, 


<t>{i)a{r) = (p{ 2 )b{r), 

<^( 2 )a(^) = (l>il)b{r), 

(E.2) 

dn(i>(i)air) = dn(i){ 2 )bir), 

dn(f>{ 2 )air) = dn(p(l)b{r), 

(E.3) 


where dn4>i'f) denotes the normal derivative of the potential functions at the branch disk, 

dn(l){r) = n-V^. (E.4) 


We further require, as a part of the boundary conditions, that the </>(*)’s be hnite at the 
boundary of the branch disk (the branch loop) and that they go to zero at infinity sufficiently 
rapidly. 

Below we reformulate the Laplace equations with given boundary conditions in terms of 
integral equations. This reformulation is known, in the case of more conventional boundary 
conditions, to be useful in many respects, including numerical computations and the study 
of existence and uniqueness properties for the solutions to the Laplace equation . We will 
use Green’s theorem, 



(r')d^r' 


for the special case where ipir') is 



- i)V<p) (r') 




(E.5) 


V^(r') 


1 


47r|r' — r 


(E.6) 
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satisfying 


V^^(r') = 6^{r' — r). 

The potential on the hrst sheet, also has a J-function source, 

- = J<53(r'-ro) 


or equivalently behaves as 

~ -f-j—/ -i’ 

^ '' 47r|r' — ro| 

when \r' — rol ~ 0. 

We define the following two important functions on the branch disk, 


(E.7) 

(E.8) 

(E.9) 


f = (Pil)a - 4’{l)b, (E.IO) 

g = 9n(/>(l)a - dn4>{l)b- (ETl) 


They are the discontinuities of the potential and the normal derivative (in the first space). 

The functions / and g may be interpreted as the dipole moment density and the charge 
density of a hctitious charge distribution placed at the disk, which, together with the point 
charge at ro, reproduces the electrostatic potential The functions / and g are the two 
unknown functions in the reformulation in terms of integral equations. 

We apply Green’s theorems to (!){{)■, choosing for M the region in obtained excluding 
the branch disk EH We then obtain 


<t>{ 2 ){r) = -x(r). 

Here we have introduced x(^) defined as 


Xir) 




1 (r' — r) 
dvr \r' — rP 


• n — 


gir') 


—-—-— 
dvr \r — r'\J 




(E.13) 


(E.14) 


where the integration variable, r’, is restricted to the branch disk. Here, the condition that 
the p(j)’s themselves be finite at the branch loop is important. Without this condition, it 
is possible to have an extra term in the above formulae originating from the branch loop. 
(Consider superposing to the 4>(^i)'s the electrostatic potential generated by a charged wire 
placed at the branch loop.) 

Eormulae (E.12)-(E.14) give an integral representation of the potential on the two sheets 
of the Riemann space, in terms of the discontinuities / and g at the branch disk. 


Substituting the integral representations back into the connecting conditions (E.2)-(E.3), 
one obtains 

Xair) + Xb(r) = -r, (E.15) 


47r|r — rol 
1 


dnXair) + dnXb{r) = 


J 


\r - rol 


(E.16) 


is also possible to avoid the use of the 5-functions. In this case one defines the manifold M excluding, in 
addition to the branch disk, the points r for both (^(i) and and ro for 4>(i)- One then requires Coulomb-like 
behaviour for ip near r and for near ro. 
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Here r is assumed to correspond to points on the branch disk. The function x also has a 
discontinuity across the disk and we denote by Xa and Xb the values of y on side a and b 


respectively. Formulae (E.14)-(E.16) are the integral equations for the functions / and g, 


which are defined on the disk. 

Eocussing on the special case in which the disk D lies in the xy-plane, one can rewrite 
the integral equations in a more explicit form as 


1 


Id v^(x - x')2 + {y- y'Y 

1 J 1 


g{x',y')dx'dy' 
J 


2 |r — ro| 


yj{x- xoY + {y- yo^ + 


(E.17) 


-1 


3e^ 


Id y y^(x - + (y - y'Y + {{x - x'Y + {y - y'Y + 

J . Zq 


= -d^ 
2 


= J- 


r — ro 


^/{x- xqY + {y - yoY + Z( 


r3- 


f{x',y')dx'dy' 


(E.18) 


They are Eredholm integral equations of the first type, the second integral equation having a 
singular kernel. We note that the hrst integral equation is equivalent to the integral equation 
for the electrostatic potential associated with a conducting disk. 


F Equations for general pp-wave matrix model 

In the main text we have given most of the formulae for the pp-wave matrix model which 
is relevant for the large J sector of the ABJM duality in the special case k = 1. In this 
appendix we give generalisations of some of the important equations to the general pp-wave 
matrix model written in terms of a mass parameter, g, and for k ^ 1. To make the use of 
the generalised formulae below more straightforward we have numbered each equation with 
the same equation number as the corresponding one in the main text, with the addition of a 
So for example equation (2.1 r) below is the generalisation of (2.1) in section]^ 
Below we consider matrices of size 




(F.l) 


We note that the angular momentum quantum number J is a multiple of k on The 

radius of R, is related to the radius of the M-theory circle, Rn, by 


^11 = f • 

k 


(E.2) 


Eor the pp-wave approximation of AdS 4 x S'^/'E^, the mass parameter g below should be set 
to the value, 

(E.3) 


Ai = 


R 


All formulae in the rest of this appendix are written in terms of the parameters AT, Ru 
and g. In order to generalise the formulae in the main text to the k ^ 1 case, it is sufficient 
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given by (F.l). 


to apply (F.2) and (F.3) to the equations below and remember that the matrix size K is 


Using the replacements X™ —)• —X"^ and U* —)• —U*, changing the sign of the mass 
parameter, —>■ —/x, and setting 2'kT = 1 brings the following equations to the form used in 
the literature on the pp-wave matrix model, e.g. [^. Note that the redefinition of /x and Y 
has the effect of exchanging instantons and anti-instantons, which can be seen from ( 2 . 12 r). 


Se = dt ti 


1 


2Rii y Dt 
Ri 


DY^ \ ^ 1 

+ 


2Rii y Dt 


DX^ 


+ (27rr)2^ ((*[X”^, X^]f + 2{i[X^, Y^]f) 


+ 


1 


(vm\2 I fJjkt 


2Rn yeJ X + 3(2.T)«„ 

_ fit’ll 




Y'‘ 


+ + 27rr— - i- ^'^ 7 ^ 23 ^ 

2 Dt 2 8 


(2.1r) 




y* = 

° 3(27rT)i?ii 


-L\ with [L\ L^] = 


(2.3i) 


r = 


/xVX 2 - 1 


(2.5i) 


Se = 


2i?ii 


dt tr 


dX” 

dt 


6(27rr)i?ii 6(27rr)i?ii 

+ {X^f + {(i[y-,y -])2 + (i[y™,y*]) 2 } 


+ ± ^y*±i(27rT)—e*^*^[y^y'=]) 

Y (st 3 Qi J 

^ X f + ii27rT)—e^^^Y^[Y\Y^]] + fermions 
V 3 3 / 


(2.8r 


T^x[y] 


+00 / \ 2 

^ V 3 ( 27 rr)i?ii J “ ^(-00)^(-00); • 


(2.10r 


iy[y] = tr l^^yV* + z(27rr)i^e*^'^y*[y^ y'^]^ . ( [2lT| r) 


^ ± ^y* ±x(27rr)—=0. 

dt 3 ^ ^ 2 L 1 J 


( 2 . 12 i) 








9 y3(27rT)i?i] 


tr(xx) = 




9 y3(27rr)/?iiy 4^ 


2 1 )x 


^^tr[(x 2 -l)l 


K,xKi\ , 


(2.15i 


where Xj = Jj/k, j = l,2,...,n, denote the (integer) size of the blocks in the case in 
which the matrices U are generators of a reducible SU(2) representation. The Kj’s satisfy 

j:,k, = k. 


= C s = 


(2.18r 
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ds 


i(27rT) 


‘iR^ ijf. 

2Cn 




(2.19i) 


pifa) 

Pi{f,9}) 



\{pif)p{9) + p{g)p{f)) ■ 




Z = 


(27rT) /6 


dvr 


dvr 

2 


K 


- I 

P 


s = 


9s 


P^ll [*^] ijk 

12 K 





d^) 


(3.3i) 


(3.6i) 


(3.7i) 


(3.8i) 
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